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We study a two-leg antiferromagnetic spin-1/2 ladder in the presence of a staggered magnetic 
field. We consider two parameter regimes: strong (weak) coupling along the legs and weak (strong) 
coupling along the rungs. In both cases, the staggered field drives the Haldane spin-liquid phase 
of the ladder towards a Gaussian quantum criticality. In a generalized spin ladder with a non- 
Haldane, spontaneously dimerized phase, the staggered magnetic field induces an Ising quantum 
critical regime. In the vicinity of the critical lines, we derive low-energy effective field theories and use 
these descriptions to determine the dynamical response functions, the staggered spin susceptibility 
and the string order parameter. 



I. INTRODUCTION 

The problem of quantum critical points (QCP) is one 
of the most important issues in the physics of strongly 
correlated electron systems, in particular in the context 
of high-T c superconductivity [Q and heavy-fermion com- 
pounds Recently this problem attracted much in- 
terest also in the context of one-dimensional (ID) quan- 
tum systems, such as ID interacting electrons, antiferro- 
magnetic spin chains and ladders, where a detailed de- 
scription of QCPs is available due to the powerful non- 
perturbative techniques based on conformal field theory, 
bosonization and integrability. It is well known that 
universal properties of ID quantum systems can be de- 
scribed on the basis of a properly chosen conformally 
invariant theory deformed by a number of perturbations 
consistent with the structure and symmetry of the un- 
derlying miscroscopic model. Quantum criticalities can 
then emerge due to the competition between two (or 
more) relevant perturbations which, when acting alone, 
would drive the system to qualitatively different strong- 
coupling massive phases that cannot be smoothly con- 
nected by a continuous path in the parameter space of 
the model. An example of a theory of this kind, display- 
ing an Ising quantum criticality, is the so-called double- 
frequency sine- Gordon model P-pl- 

In connection with ID quantum antiferromagnets, a 
plausible QCP scenario has been anticipated some time 
ago by Affleck and Haldane ||. They argued that a 
massive phase of a translationally invariant spin-chain 
Hamiltonian can be pushed towards quantum criticality 
by an external, parity-breaking perturbation. Typical 
examples of such perturbations are an explicit dimeriza- 
tion, whose role in the formation of QCPs has already 
been analyzed in several spin-chain and spin-ladder mod- 
els |7|-|l(i|], and a staggered magnetic field. While dimer- 
ization of quantum spin chains and ladders is quite re- 
alistic because it can originate, for instance, from the 



spin-phonon coupling, the case of a static magnetic field 
whose sign alternates on a microscopic scale used to be 
regarded as not achievable in experimental conditions. 

However, recently two beautiful experimental realiza- 
tions of staggered magnetic fields in quasi-lD magnetic 
insulators have been discovered. The first concerns the 
spin-1/2 antiferromagnetic chain compound Copper Ben- 
zoate |Ll[ . Due to the low crystalline symmetry, the mag- 
netic field H couples to the effective spins- 1/2 through a 
gyromagnetic tensor |l2[ 

Hmagn = E E [&P + C" 1 )" 9%] H a S& . (1) 
n a.b 

Application of a uniform magnetic field H thus induces 
a staggered field h in a direction perpendicular to H. In 
Copper Benzoate there is a second mechanism that gives 
rise to a staggered internal magnetic field. It derives 
from the staggered Dzyaloshinskii-Moriya (DM) interac- 
tion along the chain direction 

Wdm = E^ 1 )" D • x S i+!)' ( 2 ) 

3 

which, when a uniform field H is applied, induces a stag- 
gered component proportional to H x D (F|] . The pres- 
ence of a staggered field has been shown to lead to a va- 
riety of very interesting consequences [p^|-^6[ . The stag- 
gered field scenario described above is by no means spe- 
cific to Copper Benzoate. There are at least two other 
materials, Yb 4 As 3 and [PM • Cu(N0 3 ) 2 • (H 2 0) 2 ] n 
|p^| , whose properties in a magnetic field are controlled 
by the same mechanism. 

A completely different mechanism that leads to the 
generation of a staggered field has recently been dis- 
covered for spin-1 Haldane gap compounds of the type 
R 2 BaNi03 p9| , pfj| , where R is a magnetic rare earth. 
In these materials, the rare-earth ions are only weakly 
coupled to the Ni chains, but interact strongly with one 
another. They may be considered to live on a separate 
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sublattice that undergoes a Neel transition at a rather 
high temperature TV. The effect of the resulting anti- 
ferromagnetic order is to induce an effective staggered 
magnetic field along the Ni (spin-1 Heisenberg) chains 
below Xjv. 

In this paper, we study the effect of an external stag- 
gered magnetic field on the low-energy properties of the 
spin-1/2 antiferromagnetic two-chain Heisenberg ladder. 
Although at present no analogous mechanism for the gen- 
eration of a staggered field has been found for Heisenberg 
ladders compounds, in principle either of the two scenar- 
ios mentioned above are possible. One may well expect 
that a staggered field will be realized in a ladder com- 
pound before long, so that addressing this problem is 
not only of academic interest. 

The Hamiltonian of the "standard" ladder is 

.3=1,2 n 

+ J ± ]T Si," • S ^ - h E(- 1 )" 5 a,„, (3) 

n a— 1,2 n 

where J and J± are antiferromagnetic exchange coupling 
constants in the "leg" and "rung" directions, respectively. 
We employ weak-coupling ( J± <C J) and strong-coupling 
(Jj_ 3> J) approaches to show that there exists a criti- 
cal value of the staggered magnetic field, h — h c (J±, J), 
for which the system displays a Gaussian U(l) criticality 
with central charge C = 1, characterized by nonuniver- 
sal critical exponents. Both at h < h c and h > h c the 
spectrum is gapped, and the spin correlations are com- 
mensurate with the underlying lattice. This is difffcr- 
ent from the spin ladder in a uniform magnetic field |2l| l 
which induces a transition from the gapped commensu- 
rate phase (h < h c ) to a gapless incommensurate phase 
(h > h c ). Comparing the results of the weak-coupling 
and strong-coupling approaches, we find that near the 
critical point the low-energy properties of the spin lad- 
der are adequately described in terms of a XXZ spin-1/2 
chain with a J and Jx-dependent exchange anisotropy 
and an effective staggered magnetic field proportional to 
h — h c . Hence we expect that the existence of the U(l) 
QCP is a universal property of the standard Heisenberg 
spin-1/2 ladder in a staggered field. 

The critical surface h c (J±,J) separates two massive 
phases: an anisotropic Haldane spin-liquid phase at h < 
h c with coherent S z = ±1 and S z — magnon excita- 
tions having different, field dependent, mass gaps, and 
another massive phase at h > h c in which the spin ex- 
citation spectrum at q ~ n still includes coherent trans- 
verse (S z = ±1) magnons, whereas the S z — modes 
transform to an incoherent background. The transition 
is associated with softening of the S z = ±1 spin-doublet 
modes and is characterized by a divergent staggered mag- 
netic susceptibility. 

We also discuss the properties of a generalized ladder, 

Hgcn — -ffstand + V / J(Sl |7 i • Sl,„+l) (S2,n • S2,n+l) i (4) 
n 



which, apart from the on-rung interchain exchange Jj_, 
also includes a four-spin interaction V. This model is 
interesting because it can display non-Haldane, sponta- 
neously dimerized massive phases if the biquadratic in- 
teraction V is sufficiently large |2^ ]. The existence of 
such interactions in ladder compounds is supported by 
recent neutron scattering experiments p3| ]. We show 
that, in the weak coupling regime, the staggered mag- 
netic field can drive the non-Haldane phase to an Ising 
quantum criticality (with central charge C = 1/2), where 
the spontaneous dimerization vanishes and the staggered 
magnetic susceptibility is logarithmically divergent. In 
the dimerized phase (h < h c ) the spin excitation spec- 
trum is entirely incoherent, whereas at h > h c , as in 
the large-/i phase of the standard ladder, the coherent 
S z = ±1 magnons are recovered. 

The qualitative phase diagram in the (J±,h, V) space 
is shown in Fig.|l|. We note that the U(l) (Gaussian, C = 
1) and Z2 (Ising, C — 1/2) critical surfaces merge at h = 
into a critical line, which has been shown earlier |2^] 
to belong to the universality class of the SU(2)2 Wess- 
Zumino-Novikov-Witten model with central charge C = 
3/2. 
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FIG. 1. Phase diagram of the generalized spin ladder in a 
staggered magnetic field in the limit J±, V <§C J. The SU(2)2 
critical line, which at h — separates the Haldane spin-liquid 
phase from the spontaneously dimerized phase, is split by the 
field into U(l) and Z2 critical surfaces. 

The paper is organized as follows. In section 2, we 
study the phase diagram of the standard spin-ladder 
model (||) in the weak-coupling limit (J±,h <C J). Here 
we employ a field-theoretical approach which represents 
the low-energy sector of the spin ladder as an SO(3)xZ2- 
symmetric model of four noncritical 2D Ising systems 
Integrating out the fast degrees of freedom as- 
sociated with collective singlet excitations, we derive an 
effective action which describes the triplet sector with 
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anisotropy induced by the staggered field. We demon- 
strate the existence of a Gaussian criticality and show 
that, close to the critical point, the model is described in 
terms of a spin- 1/2 XXZ chain with a small anisotropy 
parameter and an effective staggered magnetic field h* cx 
h — h c . In this section we also derive projections of all 
physical fields onto the low-energy triplet sector. 

In section 3 we discuss general properties of the dy- 
namical structure factor, measured in neutron scatter- 
ing experiments, for quantum spin chains and ladders in 
a staggered magnetic field. In section 4 we determine 
the dynamical structure factor of a weakly coupled spin 
ladder for momentum transfers along the "leg" direction 
close to 7r and in both massive phases and at criticality 
In section 5 we study the induced staggered magnetiza- 
tion and show that at the U(l) transition the staggered 
susceptibility is divergent with a nonuniversal critical ex- 
ponent. 

In section 6 we consider the model (|J) in the strong- 
coupling limit (Jx)h 3* J)- Treating the exchange J as 
a perturbation and projecting the Hamiltonian onto the 
subspace of the low energy states, we arrive at an effective 
strongly anisoptropic spin- 1/2 Heisenberg chain model in 
a staggered magnetic field h* = h — J± + J/2 +■■■. As op- 
posed to the weak-coupling case, in the strong-coupling 
limit the relationship between the parameters of the orig- 
inal and the effective low-energy model is known with any 
desired degree of accuracy. In section 7 we exploit the ex- 
act integrability of the sine-Gordon model and apply the 
formfactor method to achieve a quantitative description 
of the dynamical properties of the model in the strong 
coupling limit. This section also contains a brief overview 
of the formfactor appoach. 

In section 8, we turn to the generalized ladder (Q) and 
describe the Ising transition in the non-Haldane phase, 
induced by the staggered magnetic field. In section 9 wc 
address the topological order of the generalized ladder 
model in the staggered field and analyze the field depen- 
dence of the longitudinal and transverse components of 
the nonlocal string order parameter in various parts of 
the phase diagram shown in Fig.[l| 

A comparison of the results obtained for in the weak- 
coupling and strong-coupling regimes and our final con- 
clusions are given in section 10. The paper is supple- 
mented with four Appendices which contain some tech- 
nical details used in the main text. 



II. WEAK-COUPLING LIMIT 

A. Ising-model description of the Heisenberg ladder 

In this and next sections we will be concerned with 
a weakly coupled spin ladder in a small staggered mag- 
netic field: J±,h -C J. We start our discussion with a 
brief overview of the effective field-theoretical model de- 
scribing universal properties of the spin-liquid state in 



the antiferromagnetic two- leg Heisenberg ladder |24|,|2^] . 

It is well known that a single S=l/2 Heisenberg chain 
is critical and has massless S=l/2 spinons as elemen- 
tary excitations |E5[. When a small J± is switched on, 
the spinons of the originally decoupled chains get con- 
fined to form gapped triplet and singlet excitations (for 
a review see Ref . [^6| ) . The field theory that accounts for 
the spinon-magnon transmutation in the two-leg spin-1/2 
ladder represents an 0(3)xZ2-symmetric model of four 
massive real (Majorana) fcrmions, or equivalently, four 
noncritical 2D Ising models [ p4] , p2|p^ ] . In the continuum 
limit, the Hamiltonian density 



Hm = 



es=l,2,3 



-y(»^-^ x £2)-im t ^2 



■f(&d*&-£d*£)-i™.&& 



n (5) 



describes a degenerate triplet of Majorana fields, = 
(£i>£2j £v)i ( v — L), and a singlet Majorana field ^ R L . 
The velocities Vt lS are proportional to Ja, where a is 
a short-distance cutoff of the theory. The triplet and 
singlet masses are given by 



(6) 



where A is a nonuniversal constant. The corresponding 
correlation lengths in the triplet and singlet sectors of the 
model are 



h,s ~ v ttS /\m ttS \, 



(7) 



which for a weakly coupled ladder are macroscopically 
large (lt, s 3> <*)• The last term in (|J) describes a weak 
interaction between the Majorana fermions, 



92 (t R -t L ){e R t 



(8) 



where g\ = —gi = ^TraJ±. 

For the generalized ladder fl) , the low-energy effective 
model is still of the form @, with the only difference 
that the triplet and singlet masses 



m t 



J I A - VX\ 



-3J ± X - VX' 



(9) 



(A' is a nonuniversal constant) can be varied indepen- 
dently. 

In the continuum description, the local spin densities 
of the two Heisenberg chains, Sj(x) (j = 1,2), are con- 
tributed by low-energy spin-fluctuation modes centered 
in the momentum space at q — and q — tt. Accordingly 



Sj(x) = J jR (x) + 3 jL (x) + (-l)^n,(x). 



(10) 



The chiral components of the vector currents Jj ; _R,i (i.e. 
the smooth parts of the spin densities) can be expressed 
locally in terms of the Majorana bilinears 
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K,, = J ly - J 2v = (y = R,L). (11) 

However, the most strongly fluctuating fields of the spin 
ladder, the staggered magnetizations Hj(x) and dimeriza- 
tion fields ej(x) — ► (— l) n Sj- ]n • Sj ; „+i, all with scaling 
dimension 1/2, are nonlocal with respect to £,£ 4 . These 
fields, however, admit a representation in terms of the 
order (cr) and disorder (u) op erators of the related non- 
critical Ising models [ pi] , p2|p4 | : 

n+ ~ (1/a) (mio- 2 o- 3 ^ 4 , cri^ 2 o'3A*4, cti^^^a) (12) 

n ~ (1/a) (cri^2^30-4, AilA 1 2CT3Cr 4 ) (13) 

e + ~ (l/a)/ii/i 2 /i3M4> e _ ~ (l/a)o-icr2cr 3 o-4, (14) 

where n 1 * 1 = ni ± n 2 and e^ 1 = ei ± e 2 . 

Since the the correlation lengths it ;S are large, all Ising 
systems are slightly noncritical. Whether they occur in 
the ordered or disordered phase depends on the sign of 
the corresponding mass m cx (T — T c )/T c . The crucial 
property of the standard ladder is that the signs of the 
triplet and singlet Majorana masses are always oppo- 
site. This fact, together with the known asymptotics of 
the two-point correlation functions of a noncritical Ising 
model pSfl , leads to the observation that the dynamical 
spin susceptibility of the antiferromagnetic Heisenberg 
ladder, x"(g,w), obtained by Fourier transforming the 
correlation function (n~(x,r) ■ n~(0,0)), exhibits a co- 
herent S=l single-magnon peak at u> 2 = (n — q) 2 v 2 + m| 
(with q close to n). Due to multiparticle processes, 
the dynamical spin susceptibility x"(q ~ 7r ) w ) a l so dis- 
plays an incoherent tail with a threshold at uj = 3mt- 
The singlet mode shows up only at higher energies, 
uj > 2m t + \m s \ ~ bm t . Thus, at low energies, the 
standard two-chain ladder represents a disordered spin 
liquid, similar to the Haldane phase of the spin-1 Heisen- 
berg chain with a small triplet gap. 

Let us now switch on a small staggered magnetic field 
h = hz which is assumed to be the same for the both 
chains of the ladder. The field couples to the total stag- 
gered magnetization n+, and the Hamiltonian density 
becomes 

H = Hm - (Va)cricr 2 ^3^4, (15) 

where h ~ h. Here a comment is in order. In spite of 
the already mentioned similarity between the antiferro- 
magnetic S=l/2 ladder and the spin-1 chain, it would be 
misleading to think that the role of the staggered field in 
these two cases will be similar. A weakly coupled two-leg 
S=l/2 ladder can be mapped onto a spin-1 chain by for- 
mally shifting the singlet excitation band to infinity. This 
implies the substitutions /14 — > (/X4) = 0, 04 — > (04) =/= 0, 
in which case the q = it component of the S=l spin den- 
sity is determined by the relative staggered magnetiza- 
tion of the spin ladder, n~ [see Eq.(|30|)]. So, for the S=l 
chain the magnetic term has a structure different from 



^mag = h(J,lH2<T 3 - (16) 

Since m t > 0, in the leading order the interaction ( |l6| ) 
gives rise to an effective magnetic field h ~ /i(/^i/x 2 ) ap- 
plied to the third Ising system (/1CT3). The spectrum of 
the S=l chain in a staggered field is therefore always 
massive (see Refs. 

The existence of a U(l) transition in the model ( |l5| ) 
can be foreseen as follows. Since the triplet of Ising 
copies is disordered, the magnetic term in (|l^) can be 
approximately replaced by hui cr 2 /x 4 , where h ~ h(^). 
Making a duality transformation in the fourth (singlet) 
Ising copy (m s — > — m s , [14 «-» 04), one arrives at a sys- 
tem of three disordered Ising models with the underlying 
U(l)xZ 2 symmetry, coupled by the interaction haxa^o^. 
This is an Ising-model representation of an anisotropic 
spin-1 chain close to the integrable, multicritical point 
p0[ , with a perturbation representing a parity-breaking, 
dimerization field. In the isotropic, SU(2)-symmctric 
case, this model is known |?] |l0|,^lj to exhibit a QCP 
where it becomes equivalent to the S=l/2 Heisenberg 
chain (SU(2)i WZNW universality class). A finite easy- 
plane anisotropy transforms this criticality to a Gaussian 
one. 



B. Effective action in the triplet sector 

Turning back to Eq.([l5|), we notice that the magnetic 
term contains the disorder operator /x 4 of the singlet Ising 
system which has zero expectation value and represents 
"fast" degrees of freedom of the system (although the 
ratio \m s \/m t ~ 3 may not seem large enough, it can be 
significantly increased in the generalized model (|J) with 
the "spin-liquid" condition m s m t < still preserved; see 
Eqs. (^)). We will therefore integrate the singlet mode 
out to obtain an effective action in the triplet sector. The 
existence of a Gaussian criticality will then immediately 
follow from the structure of this action. 

We write the total Euclidean action of the model as 
S = S t + S s + S st , where S t [£] and S s [^ 4 ] are the contri- 
butions of triplet and singlet sectors, respectively, and 

«,£ 4 ] = ~ / d 2 r[O h (r)+O g (r)}, 

O h = -(Va)<7 1 a 2M 3M4, O g = g 2 • £ L ) (17) 

is treated as a perturbation. Here r = (x,vt), and for 
simplicity we ignore the difference between the the triplet 
and singlet velocities. Integrating over £ 4 in the partition 
function 

Z = J D[£]D[£ 4 ] e-S&e 4 ] = const J D[|j e ~ s ^ 

yields the effective action in the triplet sector in the form 
of a cumulant expansion: 

Sefffc] = S t [£] + (S^. ~ ~ [{S 2 st ) s - (S st ) 2 s ] + ■ ■ ■ , (18) 
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where (•••)« means averaging over the free massive 
fermions £ 4 . 

The first-order correction in the expansion ( |l8| ) gives 
rise to a small renormalization of the triplet mass 

52 , / v 

m t ^m t + m s - — In , r 

Zttv \\m s \a 

The cross term proportional to g 2 h involves the correlator 
(/M( r i)£jj( r 2)££(r2))s which vanishes due to the unbro- 
ken (/i — * — /j, £ji t L —> —£r,l) symmetry of an ordered 
Ising model (see Appendix A). In the second order in g 2 , 
one obtains terms leading to renormalization of the veloc- 
ity vt and coupling constant g\ . Assuming that all these 
renormalizations are already taken into account, we are 
left with the following expression for the effective action: 



SeffK] 

^2 



4-2 /d 2 ri d 2 r 2 O t (n)O t (r 2 ) {^{vMv 2 )) t 



(19) 



where O t = a\u 2 ^. 

The fourth Ising copy is ordered (m s < 0); so the cor- 
relation function in (19) decays exponentionally at dis- 
tances r ~ L [E8fl 



(/i 4 (r)/i 4 (0)) s 



MfsY'\- r/l ,- ir/lB \ (20) 



where Ai is a nonuniversal parameter. Under the as- 
sumption that \m s \/m t 3> 1, the correlation length in the 
triplet sector, l t , is much larger than that in the singlet 
sector, l s . Therefore the product Ot(ri)Ot(r 2 ) in Eq.(pj|) 
can be treated by means of the short-distance operator 
product expansion (OPE). Using the fusion rules for two 
I sing order and disorder operators (see Appendix A, Eqs. 
([A9[),(|A10|)), we find that 



c 2 



x / d 2 r ie R e L e R e L - mi + &&) e R e L ] , (21) 

where G\ and C 2 are positive numerical constants. 

Thus, we arrive at the following effective Hamiltonian 
for the triplet degrees of freedom 



W t; eff = -■ g (Ch ■ d x £ R - i L ■ d x £ L ) - im a e R • £.1 



+ 9\\I R ll+9± {l R ll+I 2 R I 2 L ) 



(22) 



which has the same structure as the field-theoretical 
model suggested by Tsvelik Q to describe the Heisen- 
berg spin-1 chain with a biquadratic term and a single- 
ion anisotropy. The staggered magnetic field introduces 



anisotropy in the spin ladder and effectively lowers the 
SO (3) symmetry of the Majorana triplet down to SO (2) 
x Z 2 by splitting the masses, 



m d = mi 



m 3 = m t 



h 

m.« 



(23) 



and by renormalising the coupling constants of the 
marginal current-current interaction 



9\\ = 9i + C 2 v I — 



9J- = 9i 



C 2 v 



h 



(24) 



From (|23J) it follows that increasing the staggered mag- 
netic field increases the mass 7713 whereas the mass of the 
Majorana doublet, m<j, decreases and vanishes at a crit- 
ical value of the field 



h c cx \m s 



1/2 



oc jTij 11 ' 2 - 



(25) 



At this point the critical degrees of freedom are rep- 
resented by a degenerate doublet of massless Majo- 
rana fermions with a marginal current-current interac- 
tion g\\(he)- This is a typical Gaussian [U(l)] criticality 
with central charge C = 1. The vicinity of the critical 
point where the Majorana doublet (C 1 ,^ 2 ) becomes very 
soft is described by the off-critical Askin- Teller model, or 
equivalently, the quantum sine-Gordon model (SGM) for 
a scalar field $: 



Ht;eff =S ^[(d^) 2 + (d x G) 2 



m d 



cos V^nK®. (26) 



Here O is the field dual to $, the parameter 
K= 1- /2nv) + 0{g\) 



(27) 



determines the (coupling-dependent) compactification 
radius of the field and 

m d cx m t (h c - h)/h c . 

In the vicinity of the critical point, the spectral gap scales 
as the renormalized mass of the SGM (Bfl): 



M d cx - 

a \ v 



sign (m d ) 



(28) 



In section 5, we show [see Eq.( |112] )] that, in the strong- 
coupling limit (J±,h ^> J), the effective model de- 
scribing the low-energy properties of the spin ladder at 
h ~ h c represents an anisotropic (XXZ) spin-1/2 Heisen- 
berg chain with the parameter A close to 1/2 and an 
effective staggered magnetic field h* ~ h — h c . In the 
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continuum limit, this quantum lattice model transforms 
to the SGM (^) . The only difference between the weak- 
coupling and strong-coupling regimes is that in the latter 
case the parameter K is close to 3/4. 

The description of the low-energy part of the spec- 
trum of the original model ( |l5| ) in terms of the effective 
anisotropic spin-1 chain (22) holds if h <C J±. As follows 
from (|25|), this condition is satisfied at h < h c and also 
in some region above the critical field. However, if the 
field reaches values h ~ |?n, a |, the singlet mode becomes 
as important as the triplet ones, and the effective model 
( p2| ) is no longer applicable. This regime is difficult to 
tackle analytically. On the other hand, if the field h is 
further increased and occurs in the range \m s \ <C h <C J, 
the role of the interchain exchange J± becomes subdom- 
inant, and the original model reduces to two decoupled 
Heisenberg chains in a weak staggered magnetic field. In 
this case exact results are available, because in the scal- 
ing limit each such chain is described by a SGM with a 
coupling constant /3 2 = 2n (see e.g. Ref. [ p7| , chapter 
22): 



Hj = V - [(d x Qj) 2 + (c^-) 2 ] - ^ sin V2ir$j, (29) 



X(h) 



where X(h) ~ h (see |^,|l6|] for an accurate estimation of 
A(/i)). (Technically, this mapping can be achieved cither 
directly, i.e. using the rules of Abelian bosonization of 
the S=l/2 Heisenberg chain |34|||, or by establishing the 
correspondence between two pairs of the Majorana fields 
and two bosonic fields, 



the disorder operator /Z4 whose correlations are exponen- 
tially decaying at short distances, r ~ l s (see Eq.(po|)). 
Therefore one might conclude that these fields are short- 
ranged, and the spectral weight of their fluctuations is 
only nonzero in the high-energy region uj ~ \m s \. By the 
same argument, the smooth part of the relative magneti- 
zation, K, Eq. (pH) , proportional to the singlet Majorana 
field £ 4 , would also appear as a "high-energy" field. How- 
ever, this conclusion cannot be correct for the following 
reason. 

It is true that, once the high-energy singlet modes are 
integrated out, the operator n + defined in ( p2f ) has no 
projection onto the lower-energy sector. However, ( p"2] ) 
is the zeroth-order definition of this operator with re- 
spect to the staggered field which couples the high- and 
low-energy modes. In fact, apart of the always existing 
short-ranged part, the operator n + contains a strongly 
fluctuating piece, which originates from a finite admix- 
ture of low-energy states occuring already in the first 
order in h. This can be easily understood from the fact 
that reduction of the original action of the model to the 
effective one is equivalent to a unitary transformation of 
the quantum Hamiltonian of the system that projects it 
to the subspace of the low-energy states. But the same 
unitary transformation should be applied to physical op- 
erators to single out their low-energy projections. 

The low-energy projections of seemingly high-energy 
operators can be extracted from the second-order pertur- 
bative corrections to the corresponding correlation func- 
tions. Equivalently (and more formally), this can be done 
by fusing a local operator Oo(r), originally defined as a 
short-ranged field with the perturbative part of the total 
action, 



$1 ± $ 2 



O (r) - Q(r) = o (r)e 



and using formulas (A7) of Appendix A to bosonize the 
magnetic term ci 02/13/^4.) 



C. Projecting operators onto the low-energy sector 

Since the fourth (singlet) Ising system has the largest 
energy gap and stays ordered across the transition (m s < 
0), at energies u <C \m a \ the order operator 174 can be re- 
placed by its nonzero expectation value (0-4) ~ (a/Is) 1 / 8 . 
Under this substitution the relative staggered magneti- 
zation n~ and dimerization field e~, defined in (p"3|), ( p^ ) 
become projected onto the low-energy, triplet sector of 
the model described by the effective action (|2l|), (p2|): 

INT - (1/a) (cti^2M3, M1CT2M3, M1M2C3) (30) 
e- -> {a/l s ) 1/8 E-, E~ - {l/a)a x a^ (31) 

On the other hand, the total staggered magnetization 
n + and dimerization field e + are both proportional to 



-s 3t 



= 0„ (r) + l^-j j d i r 1 (0„(r)n+(r 1 )) i + 0{h% (32) 

and averaging the first-order term over the fast singlet 
modes. This term is just the low-energy projection we 
are looking for (it can be easily checked that the marginal 
part of the perturbation, given by the operator O g in 
(jl7|), yields no mapping onto the low-energy triplet sec- 
tor). This is essentially an "integrating-out" procedure 
but this time applied to the correlation functions rather 
than the action itself. 

In Appendix |b] we show that the projection of the total 
staggered magnetization onto the whole triplet sector is 
of the form: 



a 



iA r , 
iA r , 



\m s 
h 



(33) 
(34) 
(35) 



G 



where (n£) is the average staggered magnetization in- 
duced by the field (see section V), A n is a numerical 
constant and the dots stand for the high-energy parts 
of the operators. In that Appendix we also derive the 
first-order low-energy projection of the relative smooth 
magnetization of the ladder: 



K x -> -L4 



K 



K y -> -iA K 



\m s \J \l 



,1/8 / , 

at \ I l t 



\m s \J \l 



d T N-, 



K 



\m s \J \l 



1/8 



lad T .E 



(36) 
(37) 
(38) 



where Ak is another numerical constant and the fields 
N~ and E~ are defined in ©, @. 



III. DYNAMICAL STRUCTURE FACTOR IN 
THE PRESENCE OF A STAGGERED FIELD 

The scattering cross section measured in neutron scat- 
tering experiments is proportional to the dynamical 
structure factor S a °(u},Q). In this section we discuss 
some general properties of S a ^(u>,Q) in the case where 
a quantum spin chain or ladder is subject to parity- 
breaking external perturbations, such as a staggered 
magnetic field or explicit dimerization. 

Consider first the case of a single Heisenberg chain. 
The dynamical structure factor is defined as the Fourier 
transform of the spin-spin correlation function 



S^(cj,Q) 



N 



2irN 



n,m—l 



J2 / dte^-^ n - m \SZ{t)Si{0))- (39) 



Here we have set the lattice spacing <iq = 1. For a transla- 
tionally invariant, antifcrromagnetic spin chain the spin- 
spin correlation function has the following asymptotic 
structure 



<S£(t)s&(o)> 



Ff{t,n 



m) + (-l) n - m F? P (t,n-m), 



(40) 



where Fi(t,n) and F2(t,n) are slowly varying func- 
tions of n and t. According to (|Io|), in the continuum 
limit these reduce to the correlation functions of the 
smooth and staggered magnetization, ( J a (t, x) J^(0, 0)) 
and (n a (t, x)n^(0, 0)), and thus determine the dynami- 
cal structure factor in the vicinity of two different points: 
Q~0 and Q « n. 

When a staggered magnetic field is applied to a spin 
chain, the situation may seem to be different. Indeed, 
the translational symmetry of the underlying lattice is 
broken and the period of the induced magnetic structure 
is 2. The spin excitation spectrum is now defined in the 



reduced Brillouin zone — tt/2 < q < tt/2, with the points 
q = and q = ir identified. At the same time, due 
to the lowered translational symmetry, the asymptotical 
expression (^) for the spin-spin correlation function will 
contain extra oscillating pieces 

(S°(t)S&(0)) = Ff{t, n-m) + (-1)"-™^ (i, n - m) 
+(-l) n ff(t,»- to) + {-l) m Ff{t,n - m), (41) 

where F^(t,n) and F^t^n) are smooth functions that 
transform in the continuum limit to the mixed correla- 
tors (n a (t,x)jP (0,0)) and (J a (t, x)r/(0, 0)) respectively 
These correlators are nonzero in the presence of the stag- 
gered field, and the question is whether they contribute 
to the dynamical structure factor. 

The answer to this question is negative. At a formal 
level, this can be shown as follows. Note that due to the 
broken one-site translational symmetry the spin-spin cor- 
relation function should be expanded in a double Fourier 
series 



N 



(42) 



k,k' 



where k and kl vary within the paramagnetic Brillouin 
zone (—ir < fc < 7r), and the double-periodicity requires 
that Fuii 7^ if k — k' or k = k! + n. Mapping onto the 
reduced Brillouin zone yields: 



<S£(t)S&(0)>= i Yl e 



lq(n — m) 



N 



\q\<ir/2 



x[FfJt) + (-l) 



■ 7 af3 



a/3 | / -.xmTT^/S 



\7i—m TTlU 

+ (-irnT^ + (-ir^+.wi (43) 



At small q, F q ^ q , F q+7r , q+7r , F q+7! , q and F q>q +7,- are the 
Fourier transforms of the smooth functions Fi( n), FjOn), 
^3(71) and F±(n) respectively. Substituting ( [42] ) into (|39| ) 
we find that 



S a0 (co,Q) 



1 

2^ 



E 

l,V = -r 



Fq+2tI,Q+2ttV ( W ) 



(44) 



From (|4J) it follows that the structure factor is 2n- 
periodic in Q. Secondly, the r.h.s. of ( f4~I| ) does not con- 
tain off-diagonal matrix elements, Fq.q+tt and Fq+tt.q, 
implying that mixed correlators do not contribute to 
S a0 {u},Q). Therefore 



if g«o 



2tt 1 F 2 qp (w,Q-tt) if Q 



(45) 



where at small Q and to F"^(ui,Q) and F^iuj^Q) are 
Fourier transforms of the correlation functions 
(J a (t, x)J0(O,Q)) and (n a (t, x)n^{0, 0)), respectively 

The same conclusion can be reached within an equiva- 
lent but somewhat more appealing picture of diatomic 
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cells. Define a magnetic unit cell made of two sites, 
(2n — l,2n), and denote the corresponding spin opera- 
tors by S21 = Sf and S , g_ 1 = There altogether are 
four different spin-spin correlation functions 



g?(t,l-l') = (Sf'(t)S$(p)) , 
gf{t,l-l') = {Tnt) Tf(0)) 
gf(t,l-l') = {Snt) 
gf(t,l-l') = ei?(t)Sg(0)) 



(46) 



(47) 



(48) 



whose Fourier transforms 

g a (uj 7 q) = / dicxp(kjt- iq[2la }) g a (t, I) 
1=1 

have the periodicity of the reduced Brillouin zone 

g a (uj,q + Tr) =g a (uj,q) . 
It then follows from the definition (|39| ) that 

S a ^u,Q) = ^\g^ 1 (u > ,q)+gf(u,q) 

+e-^gf(cu, q) + e^gf(w, ?)]. (49) 

This expression shows that, although the spin correla- 
tion functions g a (uj,q) have the periodicity of the re- 
duced Brillouin zone, the dynamical structure factor does 
not; it rather retains the periodicity of the paramagnetic 
Brillouin zone. Thus, for the dynamical structure factor 
S a ^(oj, Q), the points Q — and Q = n/ao are inequiv- 
alent even in the presence of a staggered magnetic field. 

The functions g a (t, I) can be easily expressed in terms 
of the functions F a (t, n), defined in (f4l|). Using then fl49| ) 
one finds that 



S a f > (u,,Q) = 



N 



, dty e luJt+1Qn Ff{t,n) 

,00 N 

At e lwt+1 ^+^ n Ff(t, n). (50) 



n=l 



Using the fact that F a (t,n) are slowly varying functions 
of n, we finally arrive at the result (45) where 



N 



/ 
dt exp(icvt + ipn) Ff{t, n) . (51) 
-00 „_i 



The generalization to the case of the two-leg ladder is 
straightforward. The structure factor is defined by 



S a ?(u, q ,q ± ) 



1 



x e 



E E / d ^ 

= 1 n,m—\ 



4nN , N 

a,b=l n,m=l 
lq(n— m) — iq± (a— b) / <ja Q0 



where we have introduced a transverse momentum q± 
that can take only the two values, and n. Information 
about the low-energy part of the spin fluctuation spec- 
trum is contained in the staggered (q rs tt) and smooth 
(q « 0) parts of the structure factor for both values of 
q±. All these four cases will be considered separately. 
Below we adopt the notation 

2 2 22 
s = lo — q v , 

where q stands for a small momentum deviation cither 
from or 7r. 



IV. DYNAMICAL STRUCTURE FACTOR IN 
THE WEAK-COUPLING LIMIT 

In this section we determine the dynamical structure 
factor at low energies in the weak-couling regime. In 
subsections IV A - D it will be assumed that the stag- 
gered magnetic field is much smaller than the singlet gap 
|m s |. In this case the relevant correlation functions can 
be estimated using the effective action in the triplet sec- 
tor and the low-energy projections of the corresponding 
physical fields, discussed in the preceeding section. In 
subsection IV E we will consider another limiting case, 
|m s | -C h -C J, which will be treated on the basis of the 
model (pi). 



(52) 



A. Structure factor at qj^ = tt, q « tt 

Of primary importance is the evolution of the coherent 
triplet peak displayed by the dynamical structure factor 
of the Heisenberg ladder under the action of the gradually 
increasing staggered magnetic field. This information is 
contained in the spectral properties of the relative stag- 
gered magnetization n~ whose low-energy projection is 
given by (|30|). We therefore start our discussion with the 
case q± = tx and consider the real-space-imaginary-time 
correlation functions 

DU(r) = (a/Q^iN-^N-iO)) (a = x,y,z). (53) 

The corresponding structure factor is obtained by Fourier 
transformation, D^)(r) — > D^(q) (q = (q,(j n /v)), and 
subsequent analytical continuation to the upper complex 
w-plane (iu) n — > u> + iS). 



1. h < h c 

For fields smaller than h c the leading asymptotics of 
D^ad (r) are given by 

4)(r) = J} w (r) K (l/a 2 )(a 2 M) 1/4 

x(a 1 (r) Ai2 (r)<Ti(0)M 2 (0)), (54) 
^)(r)cx(l/a 2 )(a/^) 1 / 4 

x( M i(r> 2 (r),Ui(0W0)) (a 3 (r)a 3 (0)), (55) 
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where I3 ~ 77/7773. Here and in what follows it is assumed 
that the third Ising system (£ 3 ) decouples from the dou- 
blet (C 1 ;? 2 )- This assumption is certainly correct in the 
vicinity of the critical point and holds on a qualitative 
level everywhere in the massive phases at h =/= h c . 

Using the local bosonic representation of the product 
cti^2, given by Eq. (|A8|) of Appendix A, and properly 
rescaling the dual field O, © — > (l/\/K)Q, we find that 
the correlator in the r.h.s. of ( p4[ ) reduces to that of 
vertex operators of the dual field in the SGM (|2^) which, 
in turn, can be estimated by means of the formfactor 
bootstrap approach |37] 

/ / p — Md r j v 

-9(r) cos J ¥ e(0)) ~ Z{K) (56) 

Inspecting the r.h.s. of (|55|), we notice that the disorder 
parameters fix and [12 have nonzero expectation values 
(rnd > 0). On the other hand, since 7713 > 0, the corre- 
lator ((T3(r)<73(0)) has the asymptotic form (|20| ) with l s 
replaced by ^3. The obtained asymptotics for the correla- 
tion functions D xx \r) and Di z \r) lead to the following 
expressions for the transverse and longitudinal compo- 
nents of the dynamical structure factor: 



S+- (u, tt + q, tt) = C± {h)6 (s 2 - M 2 d ) , (57) 
S**{u>, n + q, tt) = C\\ (h)S (s 2 - m 2 3 ) , (58) 



where 



C ± (h),CJh) ~ (Im^M 2 ) 1 ^. 



(59) 



( p5| ) transform to those of vertex operators in a Gaus- 
sian model and therefore display a power-law decay. One 
easily finds that 



D£(r)oc(2) 



1/2K 



K/2 e ~r/l 3 



(60) 



After the Fourier transformation and analytic continua- 
tion we obtain: 



S + -(to,TT + q,ir) =C' A 



(A 



S**(w, w + q, tt) 



a 2 (s 2 — mi) 



6(s 2 ~m 2 ), 



(61) 



(62) 



where 



1/4 



i — Jl 

4 2 



(63) 



and 0(ir) is the Heaviside step function. Both the trans- 
verse an longitudinal susceptibilities feature incoherent 
scattering continua with thresholds at zero energy and 
lu = 7713 respectively. 



We note that the incoherent continua that contribute to 
the dynamical susceptibilities at higher energy can be 
calculated by using the results of Ref. Q (see e.g. |3(| 
for a similar calculation). 

Thus, the coherent triplet magnon peak of the isotropic 
Heisenberg spin ladder, originally located at the fre- 
quency lu — m t , is split by the field in two peaks: the 
doublet (S z = ±1) peak at u> = M d < m t and the 
S z = peak at u) = 7773 > m t . Therefore, the phase 
occuring at h < h c represents an anisotropic spin liquid 
with coherent longitudinal and transverse magnon exci- 
tations having different, field dependent mass gaps. 

Upon increasing the field, the two peaks move in oppo- 
site directions. When the critical field is approached, the 
background of multiparticle states with thresholds 3Md, 
5Md, and the doublet peak merge, and, at criticality, 
the 4-point Ising corelation functions in (Q) and ( |55| ) 
follow a power-law behavior. 



2. h = h c 

As follows from the bosonic repr e sent ation of the prod- 
ucts CTi/i 2 and (see Eqs.( A7),( A8) of Appendix A), 
at h = h c the four-point Ising correlators in (B4J) and 



3. h c < h < \m s \ 

At h > h c the mass rrid becomes negative and the dou- 
blet of Ising systems occurs in the ordered phase, whereas 
the third Ising system stays disordered (777.3 > 0). The 
behavior of D xx ' '(r) is just the same as at h < h c since 
the duality transformation associated with the sign re- 
versal of 777^, i.e. ^ix,2 ^1,2, keeps the correlator in the 
r.h.s. of d54|) unchanged. Therefore the coherent \S Z \ = 1 
magnon peak, which exists at h < h c and disappears at 
the critical point, is recovered in the h > h c phase. In 
contrast to this, the asymptotics of Z)iz (r) is changed, 
and at h > h c we find 



D(-)(r)tx JL (K|m 3 Af d 2 ) l ' ie 



-2\Md\r/v g—msr/v 



CO 



M d \r/v y/m 3 r/v' 



(64) 



As a consequence, the Haldane spin liquid loses a part of 
its coherent spectral weight at q ~ tt: the S z = magnon 
is no longer seen in the longitudinal staggered structure 
factor S zz (lu,t: + 9,tt) and is replaced by an incoherent 
continuum of states with a threshold at lu = 2\Md\ +777,3: 



S zz (lu, tt + q, tt) = C'{9[ S 2 - {2\M d \ + m 3 ) 2 



(65) 
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where 

c i'<">~(f) 



or 



1/4 



(2\M d \ 



m 3 )a 



(66) 



We will show below that, in fact, the coherent S z = 
mode still exists at h > h c but its spectral weight is 
shifted towards small momenta (q ~ 0) and is very small. 



B. Structure factor at 



q_L = 



The structure factor S aa (u, q, it) is determined by cor- 
relations of the relative magnetization 



K(x), 



1,2,3. 



(67) 



The projections of the fields K a onto the low-energy sec- 
tor are given by Eqs (p6|)-(38). 

Consider first the longitudinal structure factor. Since 
the third Ising system remains disordered across the tran- 
sition (771.3 > 0), S zz (lu, q, 7r) will be nonzero only at fre- 
quencies u > TO3. This follows from ( |38| ) and (pi]). Note 
that the operator E~ in ( |3l"| ) is related to N~ , Eq.(30), 
by a duality transformation in the doublet sector. There- 
fore the large-distance asymptotics of (E~ {r)E~ (0)) can 

be obtained from those of Di z \r), estimated in the pro- 
ceeding subsection, by interchanging the cases h < h c 
and h > h c . This leads to the following results: 
h < h r : 



S zz (u,q,TT) 

= C'{{h) ( — ) (ql s ) 2 9[s 2 - (2M d + m 3 ) 2 ]; 



(68) 



h = h c : 
S zz (u,q,n) 

h > h c 





\ ^ } 


a 2 (s 2 — m 2 ) 





S zz (w,q,w) = C ll (h) 



6(s 2 - m 2 3 ): (69) 



(ql s ) 2 S(s 2 - m 2 ); (70) 



where the prefactors C\\, Cjj, C'J are given by (59), 
(H),(H). We see that at h > h c the coherent S z = 
magnon mode is seen in the small-g part of the dynamical 
structure factor. However, its spectral weight is propor- 
tional to {ql s ) 2 and thus small. 

Consider now the transverse structure factor. Using 
the definitions (|36|), ( |37| ) and the known expression for 
the structure factor S^[oj,tt + q,ir) obtained in section 
IV A, we find that 
h ^ h c : 



S ± (u),q,ir) 

\ m., J V to, ' 



MS) 



(71) 



S (u,q,ir) 
h 

TO. 



c\ 



(72) 



where C± and C' ± are given by (£>9|), (|63|). Comparing 
the results ( |57|) and (|7l[), we see that at any nonzero h 
(except for the transition point) , the coherent transverse 
(\S Z \ = 1) magnon mode is seen both in the staggered 
(q ~ 7r) and smooth (q 0) parts of the structure factor 
S^(u), q, 7r), with the ratio of the spectral weights at q rs 
and q ~ ir being of the order of (h/m s ) 2 (M d /m s ) 2 . 



C. Structure factor at qj^=0, q « 77 

The structure factor S a ^(ui, n + q, 0) is determined by 
dynamical correlations of the total staggered magnetiza- 
tion, n + , whose low-energy projection is given by expres- 
sions (§3) ©. 



1. Longitudinal structure factor 



As follows from (p3[), S ,Z2 (w, n + q, 0) displays a broad 
continuum of states with the lowest-energy threshold 
equal to 2|Md|. If \h — h c \ <C h c , at frequencies much 
less than 7713 only the doublet modes are to be taken into 
account. In this case we have 



D^>(r)=Ai — p 



K d (r) = (e d {r)e d (0)), 



h 



K d {r), 



(73) 
(74) 



where e d = i + 1S * ne energy-density opera- 

tor in the doublet sector. A simple calculation, based on 
the assumption that the doublet fermions are free, leads 
to the result (see Appendix D): 



QmK d (q, u> + i5) = 



1 y/l 



4m 2 



2v 



(75) 



In fact, the square-root behaviour of the structure fac- 
tor near the threshold is universal: the effect of the inter- 
action between the Majorana fermions £ x and £ 2 shows 
up only in the mass and velocity renormalization and the 
interaction dependent prefactor Z d (K): 



S**(u>,n + q,0)<x±Z 2 (K) (A 
2i> V m s 



2 /, \ 2 



a 



AM 2 



[s 2 - AM 2 ) 



(76) 
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This is confirmed by the calculation done in section VI 
E, which proceeds from the bosonized correlator (|74|) 



K d (r) oc (cos V47rK$(r) cos V47rK$(0)) 



(77) 



and automatically takes into account interaction effects 
in the doublet sector within the formfactor approach 
[§7f-H to the SGM @. In that calculation, the cru- 
cial fact is that, in the range 3/4 < K < 1, the spec- 
trum of the SGM ( p^ ) consists of massive quantum soli- 
tons (s) and antisolitons (s) with the mass Md, and one 
soliton-antisoliton bou nd state (breather) with the mass 
Mi given by Eq. (127). The lat ter, however, is odd un- 
der parity (charge) conjugation ( |A16D and as a result the 
formfactor of the parity-symmetric operator cos 
between the vacuum and the breather state vanishes. 
Therefore the main contribution to the structure factor 
is due to the ss scattering continuum with a threshold at 
2|M(j|, where M d is the single-soliton mass. This explains 
the universality of the square-root behavior (|7^) of the 
structure factor S zz (ui, it + q, 0) near the threshold. 

The result ([76]) neglects the contributions of multipar- 
ticle processes with thresholds at higher energies (4|M<i|, 
6|M<z|, etc). As h — ► h c (i.e. M d — ► 0), such processes 
become as important as the two-particle ones. Exactly at 
criticality the Majorana doublet becomes massless, and 
the interaction in the doublet sector can no longer be 
ignored. In this case Kd(r) ~ (a/r) 2K , and 



C(K) 



S 2 ), (78) 



S"(u,v + q,0) oc 
where C{K) ~ [2 1 " 2 * /T 2 (K)](h/m s ) 2 (l s /a) 



2. Transverse structure factor 

The transverse structure factor S xx (ui, it + q, 0) can be 
estimated in a similar manner. The main difference is 
that now we have two Majorana fcrmions (£ and £ 3 ) with 
unequal masses (m d and TO3). Treating both fermions as 
free, simple calculations along the lines of Appendix D 
give 



S xx {oj,Tr + q, 0) oc 




(79) 



,s 2 > max {rn^m^}, 



where m± — 7713 ± m^. For different signs of the dou- 
blet mass nid ~ h c — h the frequency dependence of the 
structure factor is qualitatively different (see Fig.||). At 
h < h c {md > 0) S xx (uj,tt + q,0) follows a square-root 
increase above the threshold m + which becomes steeper 
as h — > h c . At h > h c S xx (ui,Tr + q,0) has a square- 
root singularity at the threshold m_ with an amplitude 



proportional to \m d 



1/2 



Let us now discuss the role of the so far neglected in- 
teraction between the fermions. As we have mentioned 



1.5 
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FIG. 2. Transverse structure factor at q w n, qj_ = 0. For 
more clarity, the cases h < h c (m,d > 0) and h > h c (m,d < 0) 
are shown separately. 



before, near the critical point the S z = collective modes 
described by the field £ 3 asymptotically decouple from 
the S z — ±1 (doublet) ones. Hence one can always as- 
sume that the Majorana fermion £ 3 is free and massive. 
Concerning the fermion £ 2 that belongs to the interact- 
ing doublet sector, we may use the bosoni zed e xpre ssions 
for the chiral components of £ 2 (see Eqs.( A5),(|A6| )) and 
employ a formfactor expansion in the model (^6|) to show 
that 



(&(s,r)&(0,0)) = Z{K)K X {\MA) (f) 

(£|(z,t)£!(0,0)) = Z{K)K x (\M d z\) ( - 
(&(s>r)£i(0,0)) - Z(K)K (\M d z\) 



1/2 



Z/ 
Z/ 



(80) 



Here z = r + ix/v, z — t — ix/v, and K n (x) are the 
MacDonald functions. In (§0|) only the one-particle form 
factor has been taken into account (the first correction 
involves a contribution of three-particle processes). We 
see that the expressions (|o|) have the structure of the 
2-point correlators of free massive Majorana fermions. 
The information about the interactions in the doublet 
sector, i.e. the parameter K of the SGM, is contained in 
the renormalized mass Md and the constant Z(K). This 
means that, up to this prefactor, the above result (|79| ) 
actually represent the first term of the exact formfactor 
expansion and, thus, is universal close to the threshold 
as long as h 7^ h c . 

According to (|79"|), at criticality the free-fermion ap- 
proximation leads to a step function 9(s 2 — m 2 ). How- 
ever, this result must be revisited because at Md = the 
single-fermion propagator (£ x (x, r)^ 1 (0, 0)) transforms to 
that of a spinless Tomonaga-Luttinger liquid. An esti- 
mation of S xx (lu, 7r + q, 0), quite similar to that done in 
Appendix leads to a power-law behavior 
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= C{6) 



crit 



2 2\ 28 

s z — mi 



(81) 



where 



and 



C(0) i^M ( 1 

K ' 2 1 + 2e r(l + 26) \l 3 



20 



(82) 



is the critical exponent of the single-particle density of 
states in a Tomonaga-Luttinger liquid. We see that, 
due to an "infrared catastrophy" caused by the interac- 
tions in the doublet sector, the threshold discontinuity of 
S xx (cu, 7r + q, 0) transforms to a continuous dependence. 
However, due to the smallness of 0, the power-law in- 
crease of the transverse structure factor is very steep. 



D. Structure factor at qj_ = 0, q 







At small q the structure factor S aa (uj, q, 0) is deter- 
mined by correlations of the smooth part of the total 
magnetization density. In the continuum limit, the lat- 
ter can be expressed in terms of the triplet Majorana 
fields: 



Sl,„ + S 2 ,n — > l(x), 

i a = -(i/2) e abc (£kfl+te)- 



1. Longitudinal structure factor 



(83) 



vacuum and the first breather state. This gives rise to a 
coherent delta-function peak appearing below the thresh- 
old of the two-soliton scattering continuum. Taking into 
account only the contributions of the first breather (with 
mass Mi) and the two-particle scattering states, close to 
the critical point we obtain 



S zz (uj,q,0) tx 2\- 



sin(-7r£) 1 



-47T 



Q 2 ^ 



4Mj 



1 



x exp 



3 cosh(20 o /£) + cos(tt/0 
'disinhQl -£]t)[l - cosh 2t cos(26» t/ 7r )] 



t 



sinh(2t) cosht sinh(t£) 



(85) 



where 8 = arccosh(s/2M,j) and 



A = 2cos(7r^/2)v/2sin«/2) exp 
K 



^ dt 



sini 



2 - K 



Mi = 2Af rf sin«/2) . 



(86) 



This result is valid for any 

S zz (lu, q, 0) vanishes as q - 

of the total magnetization, + S% = J dx I 3 (x), is con- 
served. In the limit K — > 1 
to ( |85| ) reduces to the form 
Md with arid. 

Exactly at the Gaussian 
tudinal structure factor is g: 
anomaly" formula (see e.g 



sign of h — h c . We note that 
because the z-component 



the two-particle contribution 
( |S4| ) as it should, if we identify 

criticality (h — h c ) the longi- 
iven by the well-known "chiral 
Ref. 



S zz {uj, q, 0) = K{a/v){qvY8{u> Z - q 2 v z ) 



(87) 



The total longitudinal current, I 3 = — IZuQi involves 
only the doublet Majorana modes which become critical 
at the transition. If the marginal interaction between 
these modes was absent, the structure factor at small q 
would display a broad continuum of states with a thresh- 
old 2\m d \ [pi: 



S**(w,?,0) 



2 2 



Am 2 d 



(84) 



As opposed to the case of the energy-density correlator 
( [77| ) where the free-fermion approximation correctly cap- 
tures its universal features, here the interaction in the 
doublet sector changes the result (|84|) dramatically. As 
already mentioned, apart from soliton/antisoliton states 
the spectrum of the SGM (§|) with | < K < 1 also 
features the first breather state. The latter is odd un- 
der charge conjugation C, wh ich inverts the sign of sine- 
Gordon field $ (see Eq.( |Al6|)). As follows from flAg| ) the 
current operator I 3 = yj K/ird x & is also odd under C 
and consequently has a nonzero formfactor between the 



2. Transverse structure factor 

Now we turn to dynamical correlations of the trans- 
verse total current, I 1 = — As before we will first 
adopt the free-fermion approximation and then discuss 
how the results are affected by the interaction between 
the doublet modes. 

A simple calculation leads to the following expression 
for S xx (u}, q, 0) (we assume that ui > 0) 

S xx (u;,q,0) 

a m 2 + q 2 (s 2 — m 2 _) + m 2 _(s 2 + q 2 )(s 2 - m\) 



2v 



l y ( s2 ~ m +)(s 2 — m ~) 



m± = 7713 ± md, s 2 > max {m^,TO^}. 

The case q = is shown in Fig.[| Apart from the 1/uj 2 
decay at ui > m 3 , the behavior of S xx (uj, q, 0) close to the 
threshold is similar to that of S xx (u,Tr+q, 0) (c.f. FigJ|). 
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FIG. 3. Transverse structure factor at q — q± — 0. 



Although at q = the transverse structure factor is 
finite for h < h c (rrid > 0), at arbitrarily small q it is 
divergent at the threshold (see Fig.[|) 



q 2 (m\ - rut) 1 /' 2 1 
m 2 + (m 2 + + q 2 y/ 2 ^ 



(89) 



where 8 — u) — 



q 2 > 0. 



For h < h c and at q ^ the curves shown in Fig.|J 
display an interesting feature. The l/y/5 drop of the 
structure factor slightly above the threshold is followed 
by an upturn which is a property of S xx (u, 0, 0). The 
maximum occurs at 



(ml 



m 2 _)q 2 



m+m_ 



At h > h c (rrid < 0) the singularity at the threshold is 
of the form: 

(to 2 -mi) 1/2 (m 2 +q 2 ) 1 

Now it survives the limit q — > and disappears only at 
the critical point. 

The free-fermion approximation is reliable as long as 
h h c . The derivation of the structure of S xx ( uj , q, ) 
at criticality is given in Appendix [5[ Using ( p9|)-(pll|) 
we find that just above the threshold id = y/ q 2 + to 2 
the structure factor S xx (uj, q, 0)| cr it only differs from the 
expression ( |sT| ) by an extra factor (1 + 2q 2 v 2 /m\) 



S xx (u,q,0) 



= C(8) 1 + 



2q 2 v 2 



■la 




(91) 



FIG. 4. Transverse structure factor at small q. 
a) md = 0.8m3; b) nid — 0.2m3. 

The discontinuity of S xx (uj, q, 0) at the threshold, ob- 
tained in the free-fermion approximation, is rounded by 
the interaction in the doublet sector. However, due 
to the smallness of 6, the transverse structure factor 
rapidly increases beyond the threshold and reaches val- 
ues of the order o f 1 exponentially close to the threshold, 
lu - y / q 2 v 2 + ml ~ oT^e^ 1 ! 26 . 



E. Summary of the structure of the magnetic 
excitation spectrum at weak coupling 

At this point it would appear to be useful to briefly 
summarize our results concerning the behavior of the 
dynamical structure factor in the weak-coupling regime 
derived above. We do this in a simple pictorial way, hav- 
ing in mind to illustrate the qualitative changes in the 
spin excitation spectrum associated with the transition 
induced by the staggered magnetic field. 

In Figs |^ and || curves denote coherent single-particle 
excitations whereas shaded areas correspond to incoher- 
ent multiparticle scattering continua. The main changes 
as we tune the staggered field through the transition 
concern S zz (uj,q 1 , n). The coherent mode visible in 
S zz (uj 7 q 7r, 7r) below the transition is relaced by an 
incoherent scattering continuum for h > h c . In the vicin- 
ity of q = exactly the opposite happens: the incoherent 
scattering continuum present below h c splits off a coher- 
ent mode for h > h c . The transverse magnon mode visi- 
ble in S + ~(u, q, 7t) is well-defined both below and above 
the transition. 
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k q 
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7i q 
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FIG. 5. Structure of low-energy magnetic excitations in 
the weak coupling limit and h < h c . 




q ± =K 



7t q 



q x =0 



o % q o % q 

FIG. 6. Structure of low-energy magnetic excitations in 
the weak coupling limit and h > h c . 



At q± = the general structure of the excitation spec- 
trum is similar for h < h c and h > h c . However, the 
precise from of the structure factor changes significantly, 
which of course has important consequences for the neu- 
tron scattering cross section. 



F. Structure factor at \m s \ <C h <C J 

This is the limit of two decoupled Heisenberg chains 
in a staggered magnetic field, described a pair of the 
/3 2 = 2tt sine- Gordon Hamiltonians (p9|). The problem 
basically reduces to the situation studied in |l3|-|l6|] . The 
elementary excitations are solitons (s) and antisolitons 
(s) with a mass gap M (which scales as M ~ ft 2 / 3 ) and 
two soliton-antisoliton "breather" bound states B\ 2 with 
gaps M and y/3M respectively. The dynamical structure 



factor has already been partially calculated by means of 
the formfactor bootstrap approach in [Q . We briefly re- 
view some impotant ingredients of this method in section 
|vi| . The leading contributions to the various components 
of the structure factor are found to be 



S+- (10, ir + q,q ± )=A ± S(s 2 - M 2 ) + 
S zz (lu, n + q, q x ) = A\\ 5{s 2 - 3M 2 ) + 

, ,2 



S+-(LJ,q,q ± )=B ± — 5( 



M 



f2\ 



S"(w,q, qi )=B\\ j^S(s 2 -M 2 ) + 



(92) 
(93) 

(94) 
(95) 



Here A±^ and B±^ are constants that can be calcu- 
lated, q±_ = 0, 7r and the dots denote higher-energy con- 
tributions of multiparticle intermediate states. Both at 
q ~ and q ~ 7r the transverse structure factor shows a 
coherent delta function peak due to one-soliton interme- 
diate states. This suggests that there may be a coherent 
transverse soliton mode throughout the whole Brillouin 
zone. The longitudinal structure factor exhibits a co- 
herent delta-function mode corresponding to the second 
breather B2 with gap \/2>M in the vicinity of q — n. For 
small momentum transfer along the chain direction the 
longitudinal structure factor exhibits a coherent delta- 
function peak which is a contribution of the first breather 
B\ with gap M. However, the spectral weight is very 
small (oc q 2 ). 

Let us now compare these results to the ones found 
above for the weak coupling limit J 3> Jj_ 3> h > h c . 



7r: The transverse stucture factor for q w n is 
similar to (^) . The longitudinal structure factor at weak 
coupling is incoherent whereas for h ^> J± it is coherent. 
This is easily understood: increasing the staggered field 
eventually splits off a bound state from the incoherent 
continuum. Finally, the structure factors for q have 
the same general form as (|94| ) and ([)5|) respectively. 

q± = 0: Here the situation at large staggered fields, 
/i> J±, is different from that found at J± 3> h > h c . 
The transverse structure factor around q = 7r and q = 
exhibits coherent modes for h ^> Jj_, Eqs. (p2]),(p4[), 
whereas it displays incoherent scattering continua for 
Jj_ S> h > h c . However these continua are singular above 
the thresholds, and one can thus easily imagine that they 
are prone to split off coherent modes once h becomes suf- 
ficiently large. An analogous situation is encountered for 
S zz (uj, q, 0) at small q. Finally, a coherent delta-function 
peak is found in S zz (uj,it + q,0) both at h 3> J± and 
Jl > h > h c . 



V. INDUCED STAGGERED MAGNETIZATION 

In this section, we estimate the magnetic field depen- 
dence of the induced staggered magnetization, M, in the 
Haldane spin-liquid phase of the spin ladder. Since at 
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h = correlations of the total staggered magnetization 
(— l) n [Sf n + S2 n ] are short-ranged, the dependence M(h) 
will be linear in the limit h — > 0: 

M(h)= Xs (0)h, (h-»0) 

X.(0) = - / d 2 r(n+(r)<(0)) ~ (96) 

On the other hand, when h is large, (h Jj_), the trans- 
verse coupling between the two S=l/2 Heisenberg chains 
of the ladder can be neglected. Each chain is then de- 
scribed at low energies in terms of the 1 — 2ir SGM 
(ph. This leads to the result 



where 



MHO x ( ' ' 



Xstag 



(h) oc h~ 2 /\ 



(97) 



To estimate the singular part of of A4(h) at \h— h c \ — > 
0, it is sufficient to know the /i-dependence of the ground 
state energy density, £(h), of the low-energy effective 
Hamiltonian. Both in the weak-coupling and strong- 
coupling limits, the latter has the structure of the SGM 
( p| ) with the amplitude of the cosine term proportional 
to \h — h c \. Using standard scaling arguments, we find 
that 

e(h)-£{h c )~-\h-h e \*^*, 



H a n = J ± S hn ■ S 2 ,„ - h(-l) n [Sl n + Sl n }. (101) 

The spectrum of the even-rung Hamiltonian consists 
of 4 levels, 



3Jj 
~4~ 



(102) 



corresponding to the triplet and singlet states 
1 



V2 



l+> = IT,T>, |o) = 
l*> = ^(IU>-|i>t»- 



I T,I) + |I,T) 



, |-> = 1 4,1), 



At ft. < J±, the singlet state \s) is the lowest-energy state. 
At h — J± the |+) level crosses with the singlet \s) and 
at h > J± becomes the lowest energy state. Focussing on 
the vicinity of the level crossing point (\h — Jj_| <C Jj_), 
we will retain two low-energy states on each even rung 



I l"|y2m — |+/2thj I ty)2m — \s)2m- 



(103) 



Projecting the even-rung Hamiltonian onto the subspace 
of these two states, i.e. imposing the constraint 
IftXlH + l Wl = l, we obtain: 



where K is given by Eq. (115). Higher-energy degrees of 
freedom which asymptotically decouple from the critical 
ones provide a finite contribution to M.(h) at h = h c . So 



M(h) = M(h, 
= m stag (h c ) + 



d£(h) 



' dh 
const I h - 



h c \*-K sign (h - h c 



The staggered susceptibility 

2(1 -if) 



Xs(h) 



\h — h r 



v 



2-K 



>0, 



(99) 



is divergent at the transition. The nonuniversal exponent 
v varies from g^/nv at J± <C J to values close to 2/5 at 
J± > J. 



VI. STRONG-COUPLING LIMIT 



In this section, we consider the model (|3j) in the strong- 
coupling limit, J±,h J. It will be assumed that the 
number of lattice sites is even, N = 2M, and periodic 
boundary conditions are imposed. 

In the th -order approximation (J = 0), the Hamilto- 
nian describes a collection of even and odd rungs, 



M 



Ho - ^ 2m 
m— 1 



M-l 



m=0 



(100) 



1 (J± 

2 I 2 



~(h-J±)U tr) 2m (tr Urn - | 4>2m(4 \lm 



(104) 



The spectrum o f th e odd-rung Hamiltonian, ffjm+i, ^ s 



obtained from (|102j) by inverting the sign of h, which 
amounts to interchanging the states |+) and |— ). Thus 
the lowest-energy states on odd rungs are 

I l"|y2m+l = |s)2m+l, I JJ-)2m+l = | _ /2m+li (105) 

and the projected odd-rung Hamiltonian reads 



PH° P---(— 



ir)2m+l(tr |2m+l - I ^)2m+l(4 bm+l f- (106) 



h-J± 



Introducing effective spin-1/2 operators T n associated 
with the n-th rung, 



T n = \(\ H)n(H\n-\ltin(ll\n\ 

t+ = i tun- \ n , t- = | H) n (it \ n , 



(107) 



we see that, in the zeroth order in J, the low-energy 
Hamiltonian describes a collection of noninteracting spins 
1/2 on a Id lattice in an effective staggered magnetic field 
h- J±: 
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PH P = const -(h- J±) 2J(-1)"T, 



The exchange interaction between the spins T„ is 
mediated by the longitudinal (J) part of the Hamilto- 
nian (^). In the two-dimensional low-energy subspaces 
of the even and odd rungs, the original spin operators 
Sj,n(j = l-i 2) reduce to 



S 



.!•» 



As a result, in the first order in J, the interaction between 
neighbouring rungs e.g. 2m and 2m + 1 is given by 



J I J 2m J 2m+l 



_J_T^S rpZ 

2 J 2m J 2m+l 



It is convenient to make a unitary transformation 



Ti± / -j \nrp± rpZ 

1 n \ 1 n ' 1 n 



rj^Z 



under which formulae (1109) transform to 



(108) coupling constants are also known to be marginally irrel - 
evan t at low energies. Finally, with the definitions ( 103 ), 
(105) of the spin- up and spin-down low-energy states, the 
effective model should be invariant under spin rotations 
around the staggered magnetic field, which is a prop- 
erty of the original model. This rules out a breakdown 
of the XY symmetry of the effective Hamiltonian, which 
would generate a finite gap. Thus, higher-order terms 
will slightly (in the order J 2 / J± ) modify the parameters 
of the effective Hamiltonian ( |112| ) and, in particular, the 
critical line (114), without affecting the criticality itself. 

(109) Small corrections to the zero-order parameter A = 1/2 
will keep the U(l) criticality far enough from the SU(2) 
critical point A = 1. 

The picture of the transition emerging in the strong- 
coupling case, J±,h^> J, is in full agreement with the 
one obtained in the weak-coupling limit, J±,h <C J (see 
section II). In both cases, the criticality is identified as 
that of an effective spin-1/2 XXZ chain. We therefore 
expect that the existence of the U(l) criticality in the 
antiferromagnetic S=l/2 two-leg ladder in a staggered 

(110) 

magnetic field is a universal property of this system. The 
vicinity of the critical field is described by the SGM ( pjjj ) 
with 



°lra 



S; 



2n 



S. 



2n 



T 
0, si 



— ( 1)™, Sf n 



>2n 



°2n 



(111) 



The low-energy effective Hamiltonian takes the form of 
an anisotropic spin-1/2 Heisenberg chain in a staggered 
magnetic field 



H cS = const + jJ2( T nT, 



L n+1 



rpyrpV 

± n ± n+l 



\rpzrpz 



where 



A = 



h* = h— Jj 



J 
2' 



(112) 



(113) 



Since |A| < 1, the model (112) has a U(l) (Gaussian) 
critical line h* = 0, i.e. 



h = J i 



^ + 0(J 2 /J±). 



(114) 



To be sure of this result, we must verify that higher- 
order terms (oc J 2 /J±), originating from virtual transi- 
tions between the low-energy and high-energy states, do 
not introduce relevant perturbations at the U(l) critical- 
ity but only renormalize the parameters of H g in (112). 
Since the original Hamiltonian (^|) is site-parity symmet- 
ric, a bond-alternating term cannot appear. Small cor- 
rections of the order J 2 /J± cannot make A greater than 
1 and drive the XXZ chain to the massive Neel phase. 
Next-nearest neighbour exchange interactions with small 



K 



1 



2(1 - iarccosA)' 



(115) 



However, the equation for the critical line, as well as 
the value of A in the effective S=l/2 XXZ chain, are 
sensitive to the strength of the interchain coupling J±. 
Being small at Ji « J, A increases with J± and tends 
to 1/2 (K — > 3/4) in the strong-coupling limit. The 

scaling law for the critical line, h c ~ J^ 2 , valid for a 
weakly coupled ladder, can be easily obtained by com- 
paring the mass gaps in the limiting cases h/ J± — > and 
Jj_/h — > 0. In the former case the mass gap is linear in 
J_l (up to logarithmic corrections). In the second case we 
have two decoupled S=l/2 Heisenberg chains in a weak 
staggered magnetic field. Since the staggered magneti- 
zation has scaling dimension 1/2, the mass gap scales as 
h~ = h 2 / 3 . The condition J± — h 2 / 3 brings us to (|2E|). 
On increasing J±, the power law h. 



J 3 / 2 gradually 



transforms to a linear dependence ( 114 ) , valid at J± J 
where the transition is governed by level crossing of the 
lowest-energy on-rung spin states. 

The strong-coupling approach can be applied to the 
generalized ladder model (Q) as well, provided that |V| <C 
J±,h. Th e eff ective low-energy Hamiltonian is again of 
the form (112), but its parameters are modified 



J — > J = J 



V 



3V 
8.7 



2(1 



8.7' 



L* , J V 

h*=h-J X+ ---. 



(116) 
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VII. PHYSICAL PROPERTIES IN THE 
STRONG-COUPLING LIMIT 

In the weak-coupling case J± <C J discussed above, 
the exact relationship between the parameters of the 
low-energy model of four massive Majorana fermions 
and those characterising the original lattice spin lad- 
der is unknown. For this reason, although the weak- 
coupling approach correctly captures the universal parts 
of all physical quantities in the vicinity of the critical 
point, the nonuniversal prefactors cannot be reliably es- 
timated. On the other hand, in the strong-coupling limit 
(J_L ^> J), the mapping onto the effective XXZ spin-1/2 
chain, Eq.( |ll2| ), is exact in the sense that all its param- 
eters can be found with any desired degree of accuracy 
in terms of the expansion in powers of J/ Jj_. Moreover, 
the projection the spin-ladder operat ors, Sf n , onto those 
of the XXZ chain, T", given by Eqs.( |lll| ), does not con- 
tain nonuniversal parameters. Therefore one can take 
advantage of the fact that the SGM (|26|), describing the 
properties of the spin ladder in a staggered magnetic field 
close to its critical value, is integrable and employ the 
formfactor bootstrap approach |37|-[4Cfl for a quantitative 
analysis of the spectral properties of the system. 

For simpli city we only consider the case V = but, in 
view of Eqs.(116), the extension to a small nonzero V is 
straightforward. 

In what follows we will use normalisations for the SGM 
that are slightly different from those used in e.g. (29), 



but are more convenient for the calculations we need to 
carry out. 

A. Scaling limit 

In the scaling limit, the XXZ spin chain with exchange 
constant J and anisotropy 

A = -costt7 2 (117) 

is described by the Gaussian model with Lagrangian 

1 



(118) 



In the strong-coupling limit we have A = ^ and thus 
, 2 



7 



The scaling limit is defined by 



(119) 



1 -7 2 

J -> oo , a -» , Ja = 2- '— = fixed. (120) 

sin tt'j 

In these conventions eto is scaled in such a way that the 
spin velocity is set to 1. It is easily restored in the fi- 
nal results by dimensional analysis. Following Lukyanov 
p0| , we will normalize the field <j> according to the short- 
distance OPE 



F - y\ 



-4 7 2 



F - y\ 



0. 



(121) 



Note that this implies that the lattice spacing must be 
taken into account explicitly when relating lattice opera- 
tors to field theory ones. For example, for the staggered 
components of the spin operators we have 



2T„ 
2T. 



± 



/ 1 \n 2-y 2 

(-1) <V 



'2A cos 



1 



(122) 



The nonuniversal constants A and F are known exactly 



F 



2(1 - 7 2 ) 2 



7 

2-2 7 a 



x exp 



dt 



sinh( 7 2 i) 



sinhi cosh(t[l — 7 2 ]) 



7 2 e- 2t 



(123) 



A = 



(2-272) 



_V^(*V). 

sinh([2 7 2 - l]t) 
sinh( 7 2 t) cosh(t[l - 7 2 ]) 



exp 



dt 

T 



-(2-> 2f 



(124) 

For 7 2 = 2/3 we obtain F w 0.5360 and A » 0.4285. 
When the staggered field term is added, the Lagrangian 
density becomes 



where = l/2 7 = and for h* < 



(125) 



2/i = - — V2Aa 



2/3 2 -l 



2T± 

n 

2T* 



F ( 1 
— exp ±1 — =( 

2 V \ 4/3 



(-l) n afW2A 



(126) 



For h* > the signs of /x and of the expression for 
need to be inverted. The spectrum of the SGM (|l25| ) at 
p 1 = 3/8 consists of soliton and antisoliton with gap Md 
and one soliton-antisoliton bound state called "breather" 
with gap 



Ah = 2M d sin(7r£/2) , ^ 



H 2 



l~f3 2 



(127) 



The soliton gap can be expressed in terms of the scale 
jj, by comparing the results of a Thermodynamic Bethe 
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Ansatz calculation with those of a perturbative calcula- 
tion valid at high energies EfJ 



n~p 2 



kT{1-(3 2 ) 



d 2 r(f/2) 



2-2/3 2 



(128) 



Combining (|128|) with (126) we can express the soliton 



gap M d in terms of the microscopic parameters of the 
lattice model 



Md 
J 



j 



4/5 



4/5 



3^(3/10) f 7rr(5/8)>/2A \ ' ^ , 
20Fr(8/lO) ^ 3V3r(3/8) J k 7 



B. Staggered Magnetization 

The staggered magnetization is given by 

((-1)" S£„) 



J + ^<(-i) n K) 



(130) 



The point to note here is that, in order to be close to 
criticality, the staggered field must be large, h fts J±, 
and consequently, the staggered magnetization is large 
as well. In the scaling limit, i.e. in the vicinity of the 
critical line h* = 0, we can use the SGM to determine 
the deviation from 1/4: 



((-l)"T n 2 )=sgn(/ 1 *) 



2f3 2 



^2A (exp(i/3</>)) . (131) 



Expectation values of vertex operators have been deter- 
mined in Ref. E2J , and in particular we have 



(exp(ip4>)) = 



(l + Qvrrq-^) 

16sm«) T{(3 2 ) 

/ r( i + |)r(i-i; 

\ 4^ 



2/3-2 



Mf 



(132) 



Combining (132) with (12C) and (124) we obtain 



<(-!)" T*) = sgn(^) 0.2897 ( L-L 



(133) 



C. Dynamical Structure factor 

Our task is now to calculate the Fourier transform of 
the retarded dynamical correlation functions in the SGM. 
This is done by going to the spectral representation and 
then utilising the integrability of the SGM to determine 



exactly the matrix elements of the specific operator un- 
der consideration between the ground state and various 
excited states. This method is known as the form fac- 
tor bootstrap approach [p8| . Let us review some of its 
relevant steps. 

In order to utilize the spectral representation, we need 
to specify a basis of eigenstates of the Hamiltonian. Such 
a basis is given by scattering states of breathers, solitons 
and antisolitons. To distinguish these, we introduce la- 
bels B,s,s. As usual, for particles with relativistic dis- 
persion it is convenient to introduce a rapidity variable 
9 to parametrize energy and momentum 

E s (9) = M d cosh 9 P s (9) = M d sinh 6 , 
E s (6) = M d cosh 9 P s {9) = M d sinh 9 , 
E B (9) = M 1 cosh 9 P B (9) = Mi sinh 9 , (134) 

where the breather gap Mi is given above. A basis of 
the scattering states can be constructed by means of the 
Zamolodchikov-Faddeev (ZF) algebra. The ZF algebra 
can be considered to be an extension of the algebra of 
creation and annihilation operators for free fermion or 
bosons to the case of interacting particles with factoriz- 
able scattering. This algebra is based on the knowledge 
of the exact spectrum and the scattering matrix |Q . For 
the SGM the ZF operators (and their hermitian conju- 
gates) satisfy the following algebra 

Z^{9 1 )Z^{9 2 ) = S e JJl(9 1 -9 2 )Z^(9 2 )Z^(9 1 ) , 

ZlW^fa) = z\{9 2 )z\,(9{)S^:k{Vi ~ 9»), 
Z^(6 1 )Zl(9 2 ) = zU9 2 )S±tU9 2 9 1 )Z<(9 1 ) 



-(2tt)£M(0i - fc). 



(135) 



Here S e ?' e ?(9) are the (factorizable) two-particle scatter- 

ing matrices and Ej = s,s, B. 

Using the ZF generators, a Fock space of states can be 
constructed as follows. The vacuum is defined by 



Z Si (0)\O) =0 



(136) 



Multiparticle states are obtained by acting with strings 
of creation operators Z\ (9) on the vacuum 

\9 n ...e 1 ) en ... ei =Zl(6 n )...Zt 1 (6 1 )\0). (137) 

In term of this basis the resolution of the identity reads 

= 1. (138) 

Inserting ( |l38|) between operators in a 2-point correlation 
function we obtain the following spectral representation 

ddi... d9 n 



^(mx^mh EE /% 



) n n\ 



x exp^Eft'^ - ejtj |(0|O(0, 0)|0„ . . . 0i> 6n ... 6l | 2 , (139) 

3=1 
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where pj and ej are given by 

Pj = M e , sinh 9j , ej = M e . cosh 9j , 

and 



(140) 



f°(9 1 . . . 0„) ei ... e „ = (0|O(0, 0)\6 n ■ ■ ■ 0i) e „... ei (141) 



are the form factors (FF). Our conventions in (14C) are 
such that M s = M s = and Mb = Mi . After carrying 
out the double Fourier transform we obtain the following 
expression for the dynamical structure factor 

i />oo />oo 

S {Lu,q) = — dx dte" i -" I (O(i,t)O t (0,0)) 

oo 



n=0 e 

X 



d6x . . . d6 n 2 

{Vl-.-V n ) ei ... en \ 

(2TT) n n\ 

5(q - ^2 M ^ sinh.0,) 5{lj - ^ Mj cosh6»j) 



(142) 



Let us now evaluate the leading contributions to ( J142) ) in 
the physically relevant cases. 



The dynamical susceptibility of the original ladder mode l 
can easily be restored by means of the relations (111) 
between the T" variables and the original spin operators 



°n,j 



s, 



,(w,7T + q,n) = S + (uj,n + q) , 



ladder 

5+ dde >,7r + g,0) = 0. 



(147) 



The nice thing about this result is that we are able 
to calculate the spectral weight in the delta- function ex- 
actly. This result is exact up to frequencies uj — 3Md, 
where additional contributions from sss states arise. 
These can, in principle, be calculated by first determin- 
ing the corresponding FFs from the relevant annihilation 
pole conditions (see e.g. [[37J for a similar calculation) and 
then by carrying out the remaining integrations over the 
rapidity variables numerically. The contribution of ( |146] ) 
to the total spectral weight at q — tt is proportional to 

4_ 

y—jj 1J an d therefore diverges as the field h* goes to 
zero. At first sight this may look strange, but the same 
is true for the gapless XXZ chain, where 



1 



(148) 



D. Transverse structure factor at q ~ n 



where n = 1 — — arccosA. 



To determine the transverse structure factor, we need 
FFs of the operator 



E. Longitudinal structure factor at q ~ n 



exp(^) 



(143) 



The lowest-lying states to which it couples are one-soliton 
states, and the first nonvanishing FF is a constant [[ll| 



|exp( —0(0,0)) |0>- = y/W) 



(144) 



In our case we have 



Zi(0) 



4 

-exp 



dtsmht sinh(t[£- 1]) 
o t sinh(i^) cosh 2 1 



r(|; 



e (i+C)t _ i 



1 



4 sinh(^) 



7 2 e" 2t 



4sinh(£i) coshi sinh([l + 



4.01 M 



2/3 



(145) 



The corresponding contribution to the dynamical struc- 
ture factor is 



F 



S+- (w, tt + q) = S(s 2 - Mj) 



0.614 



\h*\ 



8/15 



5(s 2 -M%). 



(146) 



In order to determine the longitudinal structure factor 
we need the FFs of the operator cos f3<f>. The first non- 
vanishing FF is between the vacuum and two-particle ss 
states. The FF is given by Eo| 



| cos/30(O,O)|6 



GgF(0 12 ) , 



4 smh(|6! + itt\I2Q 
dt sinh 2 (f[l - iO/it]) sinh(i[£ - 1]) 



x exp 

/(*) = 



o t 
M d ^Y{^ w ) 



sinh2t sinh£i coshi 

2/9 2 



2 r(^ 



2-2/3 2 



sinh 2 (2/3 2 £) 



2sinh(/3 2 i) sinhi cosh([l - (3 2 ]t) 



20 1 e- 2t 



(149) 



where 812 = 61—62. Carrying out the Fourier transfor- 
mation, we obtain 



A 4/3 2 |^_|2 ^ 



1 



tt s^Js 2 - AMI 



:|f(29o)f 



1.13 



J 



6/5 



7T Sy/S 



\F(26 )\ 2 . (150) 
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We note that the operator (153) reduces to the sum of the 
chiral SU(2) currents in the XXX case. Using the results 
of [|d| one can determine the leading contribution to the 
transverse structure factor at small energies by means of 
the formfactor bootstrap approach. The contribution of 
one soliton intermediate states is given by 



(J-(t,x) J+(0)) KK (M d r) 



X 2 + V 2 T 2 



K 2 (M d r), 



(154) 



where r 2 = x 2 + v 2 t 2 . Carrying out the Fourier trans- 
formation and analytically continuing we arrive at the 
following result for the dynamical structure factor 



S + (w, <7,tt) = const— o 5{s i -M£) + 



(155) 



FIG. 7. Longitudinal staggered structure factor at q± = 0. 



G. Longitudinal structure factor at q ss 



where s 2 = oj 2 — q 2 and Qq = arccosh(s/2Md). We plot 
the result in Figure 0. Just above the threshold at uj = 
\/q 2 + 4Mj the structure factor increases in a universal 
square root fashion. This is easily seen by considering 
the limit 9 < 1 of the function \F(29 )\ 2 



^ - AM 2 

\F(29 )\ 2 cx S -^L 



(151) 



Using the relations (111) we finally obtain the longitu- 
dinal structure factor of the original ladder model 



^ladder 



(w,7r + <7,0) = S zz {uj,ii + q) 



SuddcA^^ + q^ir) = 0. 



(152) 



This result is exact for frequencies below uj = 4M<j, where 
additional contributions due to intermediate states with 
2 solitons and 2 antisolitons arise. Comparing the trans- 
verse structure factor to the longitudinal one, we see that, 
in the latter case, the contribution of the ss interme- 
diate states to the total spectral weight is proportional 

2 

to (\h*\/ J)* and thus becomes small at small fields h*. 
Thus at low energies the coherent soliton modes in the 
transverse structure factor completely dominate the mag- 
netic response around the antiferromagnetic wave num- 
ber q = 7T. 



F. Transverse structure factor at q ss 

The smooth components of the transverse spin opera- 
tors in the XXZ chain are proportional to 



■exp[±^-e±04>). (153) 



The longitudinal structure factor around q — can be 
calculated analogously. The relevant Fourier component 
of the spin operator is given by 



Qq/3 
2tt 



Or 



(156) 



Here d x (j> is the topological charge density in the SGM 
and its formfactors are known |38|,[37 39 . As we have 
mentioned before, there exists a soliton-antisoliton bound 
state which gives rise to a coherent delta-function contri- 
bution to the longitudinal structure factor around q = 0. 
The contribution with the next highest threshold in en- 
ergy is due to a soliton-antisoliton scattering continuum. 
Taking only these two contributions into account we ob- 
tain 



0.0617 



J 2 ^Tm 2 



^Jq 2 + M 2 



16 q 2 y/^ 2 



AM 2 



1 



27ttJ 2 (uj 2 - q 2 )i cosh(26> /£) + cos(tt/£) 
dt sinh([l -£]£)[! - cosh2tcos(29 t/n)} 



x exp 



t 



sinh(2t) coshi sinh(i£) 



(157) 



where we again have defined #o = arccosh(s/2Af(i). The 
relation of S zz to the structure factor of the original lad- 
der model is given by (152) with ir + q replaced by q. As 
expected the structure factor vanishes like q 2 as q — ► 0. 
This behaviour is completely fixed by the Lorentz invari- 
ance of the low-energy effective theory and the fact that 
the topolog ical charge density is part of a Lorentz vector. 
The result (157) holds in the small momentum region and 
it should be possible to resolve the structure of excita- 
tions in terms of the breather bound state and the ss 
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scattering continuum by carrying out Neutron scattering 
experiments at small momentum transfer. 

Our results for the dynamical structure factor factor in 
the strong coupling limit imply the structure of low-lying 
excited states shown in Fig|| 
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q_L=o 



7i q 



71 q 



FIG. 8. Structure of low-energy magnetic excitations in 
the strong coupling limit. This picture is valid both at h < h c 
and h > h c . 



VIII. SPONTANEOUSLY DIMERIZED LADDER 

Now we turn to the generalized spin-ladder model (^) 
in which the four-spin interaction V is superimposed 
on the antiferromagnetic exchange (Jj_ > 0) across the 
rungs. If V is positive and large enough, the ladder oc- 
curs in a non-Haldane, spontaneously dimerized phase 
with a fully incoherent spectrum exhausted by pairs of 
massive dimerization kinks p2| , ff6| ]. 

The tendency towards suppression of the spin-liquid 
phase upon increasing V is al read y seen in the strong- 
coupling limit; see formulas ( |ll6|) . Within the weak- 
coupling scheme ( J±, V -C J), the transition to the spon- 
taneously dimerized phase is associated with the sign re- 
versal of the triplet mass. From Eqs. ( p"2| ) it then follows 
that, when all four Ising copies are ordered (the case 
m s ,mt < 0), the ground state is dimerized and doubly 
degenerate, with (e~) = ±e being the order parameter, 
whereas the spin excitation spectrum represents a broad 
continuum with thresholds at 2|m t | and |m t | + \m s \. 

At the transition from the spin-liquid to the sponta- 
neously dimerized phase, the Majorana triplet £ becomes 
massless (m t = 0), and the system becomes critical. The 
criticality belongs to the universality class of the level 
k = 2 SU(2) WZNW model with central charge C = 3/2 
(see Ref. p2|). When a staggered magnetic field is ap- 
plied, then, in the parameter space (m s , nit, h), the semi- 
infinite critical line m t — 0, m s < splits in the direction 
of the field h into two critical surfaces: one corresponding 



to the U(l) criticality with central charge C = 1, already 
considered in previous sections, and the other represent- 
ing the Ising criticality with C = 1/2 (see Fig.Q). The 
latter will be discussed in this section. 

The existence of an Ising QCP in the generalized lad- 
der model can be understood using an argument similar 
to the one in section II. In the spontaneously dimerized 
phase all Ising copies are ordered (mt,m s < 0). Then, 
in the leading order, the interaction term in ( jig ) can be 
replaced by /1/Z3/Z4, where h ~ h(p\o-i). If the masses 
TO3 = m t and = m s were equal, the resulting model 
would be equivalent to the double-frequency SGM in 
which an Ising QCP has already been described in much 
detail |3|Q . The existence of this transition can be easily 
visualized in the strong-coupling limit (large h). In this 
limit, the "relative" Ising degree of freedom, v = //3/i4, 
becomes effectively frozen out, while the "total" degrees 
of freedom, described, e.g., by /X3 can be tuned to critical- 
ity. This argument is still valid if the two Ising systems 
have different mass gaps, provided that they are in the 
same (in this, case, ordered) phase. 

Starting from the spontaneously dimerized phase with 
m t ,m s < 0, switching on the staggered magnetic field 
and assuming that \m s \ 3> \mt\, we can integrate the sin- 
glet mode out to arrive again at the effective model , 
with the renormalized masses still given by Eqs.([2^). The 
important difference with the previous case of the stan- 
dard ladder is that the doublet mass gap now increases 
with h while the mass TO3 decreases and vanishes at the 
same critical value as before. But since this time only one 
Majorana mode becomes massless, the criticality is of the 
Ising type. The dimerization order parameter, which is 
nonzero at h < h c , vanishes at the critical point as 

(O - {h c -h) 1/8 6(h c -h). 

The (static) staggered magnetic susceptibility XstagCO is 
constant in the zero-field limit, 

Xstag (0) ~ (l/a)(NI/KI) 1/4 (K| + KI)- 1 , (158) 

then increases with h and becomes logarithmically diver- 
gent at the transition 



Xstag(^) ~ lnfl/i - h c \/h c ). 



(159) 



Estimation of the large-distance asymptotics of the 
correlation functions is similar to what has been done 
in section 5. One only has to keep in mind that the only 
mass which changes its sign across the transition is 7713. 
Here we present the final results. 

In the dimerized phase [h < h c ) the dynamical struc- 
ture factor is entirely incoherent both in its longitudinal 
and transverse components: 



S H "U'. 7T + q. ( y 



1/4 
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S zz (uj, 7T + q, ir) cx 



At the Ising criticality 



9[s 2 - (\m d \ + |m 3 |) 2 ] 
y/s 2 - (\m d \ + \m 3 \) 2 ' 

%y g[a a -4m3] 
W3 / ^/s 2 - 4m 2 , 



(160) 
(161) 



5* (w,7r + g,7r) oc ( — 



„2 \ 4 



?[a; 2 — q 2 v 2 — m 2 ,] 
\s 2 — m 2 ,] 3 ^ 4 



S zz (to, 7T + q, 7r) cx 



W 2 L 



[s 2 - 4m 2 ] 4 



(162) 



(163) 



In the region h > h c , the longitudinal staggered spin 
fluctuations remain incoherent: 



S zz (u>, 7T + q, 7r) cx 



^3^2 + ^ 3 ) 2 



2\ 4 



x 9[s 2 - (2|m d | +m 3 ) 2 



(164) 



whereas the transverse part of the dynamical structure 
factor displays a coherent 5-function peak, 



S + (u), n + q, 7r) cx 



S (s 2 - m 2 ) , (165) 



describing a massive magnon with the spin projection 
\S Z \ = 1. As in the standard ladder, at h > h c the 
dynamics of e_ also displays a coherent mode with the 
mass 1 777,3 1. 

Thus, characterization of the massive phase occuring 
at h > h c in the spontaneously dimerized ladder coincides 
with that for the standard ladder. This phase occupies 
the region separated by the U(l) and Ising critical sur- 
faces on Fig.|l|. 

IX. STRING ORDER PARAMETER 

Den Nijs and Rommelse [^J and Girvin and Arovas 
p8[ have shown that the Haldane gapped phase of the 
spin-1 chain is characterized by a nonlocal topological 
string order parameter, 



m— 1 
j=n+l 



(O a )= lim (S£exp(i7r X S ?) S nr)> (a = x,y,z), 

\n—m\ — >oo 



whose nonzero value is associated with the breakdown of 
a hidden Z2 x Z2 symmetry p9| . For a weakly coupled 
SU(2)-symmetric spin-1/2 Heisenberg ladder, the string 
order parameter, defined as 



°n,m = II (- 4S ?i^) = 6X P 



2:P 



(166) 



was discussed in Refs. p4,p2| (see also p7[|). The de- 
scription of the low-energy degrees of freedom of the spin 
ladder in terms of the Ising variables i s esp ecially effi- 
cient in this case because the operator ( |166| ) acquires a 
simple local form in terms of the Ising operators er Q , /i a 
(a = 1,2,3), and the hidden Z2 x Z2 symmetry becomes 
manifest. In the continuum limit [p0[ 



lim (O a (x,y)) = (O a ) ~ (073 £r 7 ) + (npn^) , 

\x— y I — * oc 



(167) 



so that the string order parameter reveals the SU(2) sym- 
metry and is nonzero both in the Haldane and dimerized 
phases just because the degenerate triplet of the Ising 
systems is either disordered ((<r a ) = 0, (/x a ) 7^ 0) or or- 
dered {(a a ) 0, ( M a) = 0). 

As we have seen, the staggered magnetic field removes 
the SO (3) degeneracy of the triplet modes, and in the 
massive phase located between the two critical surfaces 
shown in Fig|l] (the case h > h c ), the signs of the masses 
mi = 772 2 and 7773 are opposite. As a result, the depen- 
dence of the longitudinal and transverse components of 
the string order parameter on h becomes qualitatively 
different. 



A. U(l) transition in the Haldane phase 

Let us start with the longitudinal string order param- 
eter, O z . In the region of small fields, < h <C h c , one 
can adopt the picture of three independent triplet Ising 
copies as the zero-order approximation and take into ac- 
count the effect of the staggered field as a small pertur- 
bation. According to ( p3|) and (pij), the latter leads to 
splitting of the triplet masses and renormalization of the 
coupling constants, both being of the order of h 2 . As a 
result 



(O z )(h)c(^)( m ) 

= (a)(o) 



Cx {I 



h 

m. 



0(¥) 



(168) 



where (O z )(0) ~ (rrita/v) 1 ^ 2 . In vicinity of the critical 
point, \h — h c \ <C h c , the Ising doublet (1,2) becomes very 
soft and asymptotically decouples from the rest of the 
spectrum, being described in terms of the SGM (p6|). In 
this case the operator O z can be bosonized (see Appe ndix 
A). However, since if ^ 1, in formulas ( |A7| ) and (AS) 
for the products of Ising operators and cxi<72, one 

should rescale the field <&: <& — > \/K$. So 
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(O z ) ~ (smV^K®) 2 + (cosV^$) 2 . (169) 
Since ~ h c — h, at h < h c (sin \J i:K<&) — 0, and 



(cos -vArfC $) ' 



M. 



AT/2 



(/lc-/l) = 



(170) 



At h > h c (cos V7rX$) = 0, and 



\M 



K/2 



(h-h c ) — . (171) 



Formulas (17C), ( 171 ) determine the power law according 
to which the string order parameter (O z ) vanishes at the 
transition. 

Upon increasing the field in the region h > h c (O z ) 
keeps growing. In the limit of strong fields, h 3> h c , 
the system represents two identical copies of the Heisen- 
berg S=l/2 chain in a staggered magnetic field, each 
of them being represented in the continuum limit by a 
f3 2 = 2tt SGM with the nonlinear term proportional to 
hcos y/2w$> a (a — 1,2). Clearly, the field $ is a symmet- 
ric linear combination of <I>i and $2: 



$ = 



$1 + $2 



Therefore, in this limit 



<"M/f$i) 2 (cosy!$2> 
(sin ^|# 1 ) a (sin y~$2> : 
(co S y|$i) 2 <sinW|$ 2 ) 
(sin J ~$x) 2 (cos J |$ 2 ) 



(172) 



The minima of the potentials /i cos v 27r<I> a are ( < &a) m = 
V2n(m + 1/2) at h > and ($ a ) m = V2irm at h < 0, 
wh ere m — 0,±1,±2, .... Therefore, the last two terms 
in ( |172| ) vanish, and for any sign of h the z-component 
of the string order parameter grows as 



(o.) ~ H 1/3 . 



(173) 



Consider now the transverse components of the string 
order parameter, (O x ) = (O y ). At h <C h c the behaviour 
of (O x ) is similar to that for (O z ): 

= (a)(o) 



4 U 



9i + c aluhl 

7T TTV Vlt 



h 

m., 



hi 



0(h 4 



(174) 



Near the critical field, where the doublet of the Ising 
systems decouples from the third Ising component of the 
triplet, one finds that 



(O x )(h) 



M d a 



1/4 



9(h c - h) 

(h c - h)^^e(h c - h). 



(175) 



Thus, due to the fact that at h > h c the Ising doublet 
becomes ordered ((^1) = {^2) = 0), in the large-/i mas- 
sive phase the transverse components of the string order 
parameter vanish. 
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FIG. 9. String order parameters as functions of the stag- 
gered magnetic field across (a) the U(l) transition in the Hal- 
dane phase and (b) the Ising transition in the dimerized phase. 



B. Ising transition in the dimerized phase 

In the spontaneously dimerized, non-Haldane phase 
the Ising doublet remains massive at any h. The mass 
gap |m<j| is an increasing function of h; so is the longitu- 
dinal string order parameter. At h <C h c 



(O z )(h) ~ (oiXaa) 
= (C)(0) 



2 \a 



— +0(h 4 



(176) 
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(O z ) grows monotonically with h and at large h (h^> h c ) 
crosses over to the {hi 1 / 3 be haviour. At h <C h c (O x ) is 
still given by formula (174), whereas close to the Ising 
criticality 

(O x ) ~ (h c - h)^ 4 d{h c 



h). 



(177) 



The dependence of the string order parameters on the 
staggered field is schematically shown shown in Figs.^. 



X. CONCLUSIONS 

In this paper, we have analyzed the properties of the 
two- leg antiferromagnetic spin- 1/2 ladder in a staggered 
magnetic field. We have considered the spin-liquid phase 
of the standard ladder and the spontaneously dimerized 
phase of a generalized spin-ladder model. We have shown 
that in the former case the staggered field drives the sys- 
tem towards a Gaussian criticality, whereas in the lat- 
ter case it induces an Ising transition. These two criti- 
calities are associated with a softening of the transverse 
(S z — ±1) and longitudinal (S z — 0) collective modes, 
respectively, and are characterized by power-law and log- 
arithmic divergencies of the staggered magnetic suscep- 
tibility. 

By comparing our results for weakly (Jj_/ J -C 1) and 
strongly coupled (J±/J ^> 1) ladders in a staggered field, 
we can identify certain universal features of the transi- 
tion. The very existence of the Gaussian criticality is uni- 
versal. At criticality the low-energy degrees of freedom 
can be decribed in terms of spin-1/2 spinons similar to 
those found in the anisotropic spin-1/2 Heisenberg XXZ 
chain. This is interesting and shows that the staggered 
field h leads to a destabilisation of the magnons that form 
the low-lying part of the spectrum in the absence of h and 
eventually "deconfines" them into pairs of spinons. 

Close to this criticality the low-energy part of the spec- 
trum involves a doublet of interacting transverse modes. 
These are described in terms of an effective spin-1/2 XXZ 
chain in a weak staggered field, which vanishes at h = h c . 
The transverse modes are identified with quantum soli- 
tons (with gap Md) of an underlying SGM. These solitons 
together with soliton-antisoliton bound states determine 
the behaviour of the dynamical structure factor in both 
massive phases (h < h c and h > h c ). The soliton reveals 
itself as a coherent delta-function peak in the transverse 
staggered structure factor S ± (cu, n + q,q± = ir). 

On the other hand, the q± = part of the spin excita- 
tion spectrum is different in the weak and strong coupling 
cases. The same hold true for the q± = ir part of the lon- 
gitudinal structure factor: at weak coupling the S z — 
mode is still seen in the low-energy part of the spectrum, 
whereas it is pushed to very high energies in the strong 
coupling regime. 

In the absence of a staggered field the generalised lad- 
der is spontaneously dimerized and its elementary exci- 
tations can be understood in terms of topological dimcr- 
ization kinks p2[ . When a staggered field is applied the 



dimerization diminishes but the qualitative picture re- 
mains unchanged until the field reaches its critical value 
h Cl where we find an Ising criticality. Here once again 
we may think of the elementary excitations in terms of 
pairs of spinons. The physical properties for h > h c are 
the same as in the h > h c phase of the standard lad- 
der discussed above (see also Fig.0). In particular, the 
coherent S z = ±1 modes are recovered. This identifi- 
cation is confirmed by the behaviour of the longitudinal 
and transverse string order parameters. 

As we have seen, the physics of the spin-1/2 ladder in 
a staggered field is very rich. We think that it would be 
very interesting to explore it experimentally. An open 
question is to analyze the crossover region between weak 
and strong coupling regimes by e.g. numerical methods. 
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APPENDIX A: SOME FACTS ABOUT 2D ISING 
MODEL 

In this Appendix we briefly summarize those facts 
about the 2D Ising model which are used in the main 
part of the paper and other Appendices. 

Close to criticality (\T — T c \ <C T c ), the scaling prop- 
erties of the 2D Ising model are described by a Lorentz- 
invariant (l+l)-dimcnsional quantum model of a mas- 
sive real (Majorana) fermion p^|3^]. The correspond- 
ing 2D Euclidean action (written in complex notations: 
z = T + ix, z — t — ix, d — d/dz, 8 = d/dz; v = 1) reads 
ill: 



S = d 2 z fadh + Z R dS R + vmZUn) 



(Al) 



Here £l and £r are the holomorphic (left) and antiholo- 
morphic (right) components of the fermionic field. The 
magnitude of the mass, 



T-T c 



(a being a short-distance cutoff) determines the correla- 
tion length in the Ising model, l c ~ u/|m| 3> a, which 
diverges at criticality (m = 0), and the sign of the mass 
indicates whether the system is ordered (m < 0) or dis- 
ordered (m > 0). The set of strongly fluctuating fields of 
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the Ising model (at criticality these are known as primary 
fields of the conformal field theory with central charge 
C = 1/2) includes the fermion field (£l,£r), the mass 
bilinear (or energy density) e — i£r£li an d order and 
disorder fields, a and \x. The latter two fields are non- 
local with respect to each other and also with respect to 
the Majorana field. 

Two identical noninteracting copies of the 2D Ising 
models are described by a pair of Majorana fermions, £ x 
and £ 2 , which can be combined into a complex (Dirac) 
massive fermion, ip = (£* + i£ 2 )/\/2. The latter can be 
bosonized ||^] . If the two Ising copies are slightly noncrit- 
ical, the resulting bosonic theory represents a quantum 
sine-Gordon model at the decoupling (or Luther-Emery) 
point f3 2 — Air: 



S SG = [ d 2 z i(3 M $) 2 - — cos 
J \_2 na 



(A2) 



Below we give a list of bosonization rules for two Ising 
copies |3^,g7],|| which are used in the main text: 

• The chiral components of the U(l) current: 



=d(j) L {z), 



me) = ie R (mm = -^bm*)- (A3) 

\/7T 



The total energy density: 

£1(2,2) +£2(2,2) = — cos VAtt<&(z,z). (A4) 

ITCH 



• The fermionic fields: 

- (7ra)- 1 / 2 exp[-i^ i (2)] ) (A5) 
+ - (^a)- 1 / 2 exp[iv/4^ fl (z)]. (A6) 

• Mixed products of the order and disorderd opera- 
tors (a more accurate definition of these products 
includes Klein factors &): 

a 1 (j2 ~ sin y/TT$, fiiH2 ~ cos v / tt < I ) , (A7) 
C1M2 ~ cos-x/tt©, Mi ct 2 ~ sm^f^^Q. (A8) 

In the above formulas, 

$(2,2) = <fo(z) + <Mz) 

in the scalar field of the underlying SGM, and 

6(2,2) = cj) L (z) - cf) R (z). 

is its dual counterpart. 

Using these bosonization rules, one can easily recover 
all OPEs for a single Ising model. In particular, fusing 



the products of Ising operators in (A7) one derives the 
OPEs 



a(z, z)a(w, w) 



1 



V2 \\z - w\ 

/Lt(2, Z)fi(w, w) 
1 



1/4 



7T 2 



w\e(w, w)} 



V2\\z 



w 



1/4 



[1 — tt\z — w\e(w, w)} 



(A9) 



(A10) 



used in section II B. In the same way, from representation 
(|ASD one derives two more OPEs: 



a(z, z)fi(w, w) 

2 (|z-w|/a) 1/4 
li{z, z)a(w, w) 

Fk_ j*(z - w)^ 2 Cl(w) + 7(2 - w) 1 / 2 ^^) 
2 (|z - iu|/a) 1/4 

where 7 = e l7r / 4 . 

Since the SGM ( |A^ ) occurs in a topologically ordered, 
massive phase, in the ground state the field <& is locked in 
one of the infinitely degenerate minima of the potential 
U = — to cos yJ\~TT<&: 



(AH) 



(A12) 



($); i ac = if m > 0, 

($)™ c = ^(n+1/2), if m<0. 



From (A7) it then follows that 



(0-1,2) = 0, (^1,2) ^0, if to > 0, 
(ai, 2 ) ^ 0, (mi >2 ) = 0, if to < 0. 



(A13) 



(A14) 



Quantum solitons of the SGM ( [A2] ) are associated with 
the vacuum-vacuum transitions, <I> — ► $ ± ^/tt, O — > O, 
which correspond to the following Z2 x Z2 transformations 
of the fermionic fields and Ising operators: 



Si — * Si! Si? 
(Tl — > =F<71, CT 2 - 



^-^ (a = 1,2), 
±fj 2 , Mi ±Mi, 



TM2 (A15) 



For any sign of to, this symmetry is spontaneously broken 
in one Ising copy and preserved in the other, leadin g to 
the conclusion that quantum solitons of the model ( |A2| ) 
describe kinks of a single ordered (disordered) Ising sys- 
tem that connect opposite values of the order (disorder) 
parameters. 

Parity (or charge conjugation) transformations 



to > 



f 2 

Sfl,i 



4 £fl,i! 

e->-e, 

fl<T2 — ► — 0"lO"2, MlA*2 
CT1/X2 -* CTi/i 2 , /il(T 2 — 



MlM2, 
"MlO-2 



(A16) 



and 



25 



0102 — > (7i(T2, MlA*2 _ * —^1^2, 
(TlM2 ~> ~CTlM2, M1CT2 -> A* 1^2 



(A17) 



keep invariant vacuum expectation values of the order 
(disorder) parameters at m < (m > 0) and therefore 
serve as a tool to conclude whether a given correlation 
function in a broken-symmetry phase is nonzero. For 
example, consider the correlation function 



K (r) = (M(r)e(O)) 



(A18) 



for a single ordered Ising model. For two identical and 
decoupled Ising models this correlator can be squared: 

K 2 (r) = ( Ml (r)/x 2 (r)ex(0) £2 (0)) = K 2 {v). (A19) 



Under transformations (A17) the product £162 stays in- 
tact but changes its sign. Therefore 



K 2 (r)=K 2 (r)=0. 



This fact has been used in section IIB. 



(A20) 



APPENDIX B: FIELD INDUCED ADMIXTURE 
BETWEEN THE SINGLET AND TRIPLET 
MODES 

In this Appendix we consider two, apparently "high- 
energy" , operators: the total staggered magnetization, 
n + , and the smooth part of the relative magnetization, 
K, defined in ( |l2|) and (|ll|), and find their projections 
onto the low-energy, triplet sector of the model. 



1. Projecting n + 

As follows from the comparison of (|3^) with formula 
(|l9|), the low-energy projection of the operator O = nf 
has actually been found in section IIB, and the result is 
contained in the second-order correction to the effective 
action Thus we arrive at Eq.(|33|). 

Consider now the operator Oq = nj. Treating all Ising 
systems as decoupled, we have: 



(ni(r)ni (rt)) s = a 2 [/ii(r)oi(ri) 
x [o- 3 (r)// 3 (ri)](//4(r)/i4(ri)) s 



[CT 2 (ri)cr 2 (ri 



The correlator (/X4(r)/i4(ri)) s is short-ranged (see (p0[) ) . 
Hence the products of operators in the square brackets, 
all of them defined in the triplet sector, are subject to 
fusion. Using the fusion rules (AE), (All) and ( A12| ) 
and integrating over the relative coordinate p = r — ri , 
we arrive at the result (|34|). In the same way one obtains 
the low-energy projection for n+ given by j35|). 



2. Projecting K 

Here we derive the low-energy projection of the oper- 
ator 

K = i(^ + ^0- 

Consider first the case Oq = K x . The correlator in (|3 
reads: 

(^(r)n+(ri)) s = (i/a)a 2 (r 1 )fi 3 {r 1 ) 

x £ ^(r)a 1 (r 1 )(^(r) M4 (r 1 )) s (Bl) 

v=R,L 

Keeping in mind that m s < 0, first we note that the corre- 
lators '■{£^(j)^4 : (ri)} s are invariant under charge conjuga- 
tion ( A17) and, therefore, are nonzero. These correlators 
are chiral but otherwise short-ranged, decaying exponen- 
tially at |r— ri| ~ l s . Since they serve as integral kernels, 
a qualitatively correct estimation of the integral in ( |32] ) 
can be obtained if one treats the product 

*i(*i,*iM*i,*i) (B2) 

by OPE and then confines the integration region to the 
interval < |r — ri| < l s . To proceed further, one can 
bosonize two local products, £*£ 4 and <7i/i4, and then 
fuse these fields as those belonging to a (critical) Gaus- 
sian model. Since we are looking for a short-distance 
OPE, the relative sign of the Majorana m asse s mi and 
TO4 is unimportant, and bosonization rules (A7), ( |A§| ) are 
perfectly applicable. We obtain: 



1 



[dcj) L (z) - d(j>n{z)] cos y/nQ (z 1} zi) 

1 1 



1 

4-7T 



Z — Z\ z 

1 $le(z — z\ 



Zl 



sin y/Or&fa, zi) 



2tt \z 
and therefore 



fii(zi,zi)a 4 (z,zi), 



(X x (r)n+(ri)) s 



Tl 



2?r Ir-rJ 



(B3) 



(B4) 



where N~ ~ a^ 1 ^iiG2H 3 {(Ja) is the y-component of rel- 
ative staggered magnetization averaged over the high- 
energy singlet modes. According to the definition ( p2| ) 



where 



r — ri 



(PQ)Pl) = (vt,x). 



Expanding in p, the lowest-order projection of K x onto 
the triplet sector is found to be given by formula (|36|). 



2G 



Quite similarly one arrives at formula ( |37j ) for the pro- 
jection of K y . 

In the cas e Oq — K z , we follow the same procedure 
and use (A7) to bosonize the product ^3^(4. The corre- 
sponding OPE reads: 

~ \p(j}{z) ~ d(f>(z)] COSy/lT^(z 1 ,Z 1 ) 



1 



1 



1 



4ir \ z — z\ z — z\ 
i 5m(z — z±) 



sin t/tt$(zi, Z\) 



2n \z — Z\\ 



a a (zt, z\)ai{z, zi), 



and this eventually leads to formula (pq) for the low- 
energy projection of K z . 



APPENDIX C: ENERGY-DENSITY 
CORRELATIONS IN THE DOUBLET SECTOR 

Here we estimate the Fourier transform of the 2-point 
energy-density correlation function in the doublet sector, 



K d (q,LUm) 



J dxdr K d (x,r)e 



— l(qx—€r) 



(CI) 



K d (x,T) = {T T e d (x,T)£ d (0,0)), e d 



a=l,2 



and find its analytical continuation (iw n — ► uj + i<5) that 
determines the structure factor S zz (oj, ir + q, 0). We will 
consider the case 7^ and assume that the two 
massive Majorana fermions are decoupled. In this case 
K d (r) = 2 detG(r), where G(r) is the real-space 2x2 
Green's function matrix for a free massive fermion in the 
Nambu representation. We have (see Fig. fil]) 



Zir Zir 



LL\ 



- G RL (k + ,E + )G LR (-k-,-e-) ], 



(C2) 



where 



G(fc,e) = 



G RR (k,e), G RL (k,e) 
G LR [k,e), G LL (k,e) 

ie + kvr 3 + m d T 2 
e 2 + k 2 v 2 + m 2 d ' 



(C3) 



and k± — k± q/2, e± — e ± w„/2. 

Since the mass bilinear is a Lorentz-invariant ob- 
ject, K d (q,uj n ) depends only on uj 2 + q 2 . This makes it 
sufficient to estimate K d (0,uJ n ) in which case the calcu- 
lations become especially simple. Integrating over e in 
@ yields: 



K d (0,u n ) 



dk 1 

27r iw„i? fc 



iWnE'fe - 2m 2 , iio n E k 



2m 2 



iuj — 2Ek 



2E k 



(C4) 



(x,x) 




(0,0) 



FIG. 10. Two-fermion bubble. 



where Ek — \J k 2 v 2 + m\. Analytically continuing this 
expression, then taking the imaginary part and finally 
replacing to 2 by s 2 = oj 2 — q 2 v 2 , we get (ui > 0): 



$smK d (q, lu + id) 



1 y/ s 2 — 4m 



2v s 

This result has been used in Eq.([76|). 



(C5) 



APPENDIX D: TRANSVERSE STRUCTURE 
FACTOR AT SMALL MOMENTUM 

In this Appendix we provide some technical details 
concerning the small-g structure factor S xx (oj, q, 0) at the 
critical point. We start from the real-space representa- 
tion of the Matsubara polarization operator X xx (x, r) 



X**(x,t) = (TtI^x^t)! 1 (0,0)) 
= Tr G(t,x)Q(t,x) 



(Dl) 



Here G and Q are is the Green's function matrices for the 
massive and massless Majorana fields £ 3 and £ 2 , respec- 
tively. Due to the marginal interaction in the doublet 
(C 1 ;^ 2 ) sector, Q has the structure of the single-particle 
Green's function of a spinless Tomonaga-Luttinger liq- 
uid with the interaction constant K. Using the explicit 
expressions for these two Green's functions, 



G RR (z,z) 
Gll(z, z) 



m 3 z 



2ir V z 

777,3 jl 



2tt V z 



Ki(m 3 \z\), 
Ki(m 3 \z\), 



(D2) 



Grr(z,z) = -_(-) ( — 

Gll(z,z) = - 
G R .l(z,z) = 0, 
we obtain: 
X(z,z) 
m 3 a 



1 /z\ 1 / 2 I n 

2na \z/ V l z L 
2na \z) \ \z\ 



(D3) 



(27T«) 2 



±{K+1/K) 



V z z/ 
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Introducing two-dimensional real vectors, q = (uj n ,qv) 
and r = (r,x/v) (u>„ being the Matsubara frequency), we 
pass to the Fourier transform X xx (q) and, after angular 
integration, obtain: 



X xx (q) = --(m 3 a) 



20 



: / dx x 
'o 



-26 



Ki(x)J 2 X 

1 m 3 



where J 2 {x) is the Bessel function, and 



29 =2 



K+ K 



(D5) 



(D6) 



Since th e ex pressions for the Tomonaga-Luttinger prop- 
agators (D3) are asymptotic (i.e. valid at \z\ > a) the 
lower integral cutoff must be finite, ~ m^a. On the other 
hand, the integral in ( |D5| ) is convergent at a; — > 0, if 
6 < 1. As follows from (D6), this condition is satis- 
fied not only in the weak-coupling case, where K is very 
close to 1, but also in the strong-coupling regime where 
K = 3/4 (in the latter case ps 2 • 1CT 2 ). This justifies 
the replacement of the lower integral limit by 0, in which 
case the result can be expressed in terms of a hypergeo- 
metric function: 



X xx (q) 



1 r(2 - 0)r(l - 0) 



2 1 + M r(3) 



(m 3 a 



,20 



? 1 

< - <r 



F 1-0, 2- 0; 3;--^ 



(D7) 



To single out the leading singularity at the threshold in 
the interacting case, we use the transformation formula 
f32l, 



F(a, b; c; z) 
T(c)T(c-a 



b) 



r(c- a)T(c-b) 

_ a _ b r(c)T(a + b-c) 



F(a,b;a + b- c+ 1; 1 - z) 



+ (1-*) 
xF(c — a 



T(a)T(b) 
b:c — a — 6+1:1 



and formally treat 1—z — l+q 2 /m| as a small parameter. 
Then in the leading order 



X xx (q) 



i r(-20) 

„2 N 2e 



(m 3 a) 2e 



^7, 



1 



rn , 



1 + 1 



(D9) 



Performing analytical continuation (iw„ 



iS) 



{tul + q 2 + mj) 



20 



rn 



2 \ 29 [9(m 2 3 - s 2 ) + 6{s 2 - ml) cos 2tt0] 



,2|20 



9{s 2 — ml) sin27T0 signw 



(D10) 



and using the relation 

r(-20) sin27T0 = -, 

v ' r(i + 20)' 

we arrive at the expression (^l|) of section IV D 2 
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